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figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 


Xll 


Fig.  15b. 


Fig.  15c. 


Fig.  15d. 


Fig.  15e. 


Fig.  15f. 


Fig.  15g. 


Fig.  15h. 


Fig.  15i. 


Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  113232,  at 
8.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  16a. 


Fig.  16b. 


Fig.  16c. 


Fig.  16d. 


Fig.  16e. 


Fig.  16f. 


Fig.  16g. 


Fig.  16h. 


Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  op  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  axe  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  17a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1123013, 
at  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right- 
hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 

Fig.  17b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line.  Cross-sectional  plot  given  in 
Figure  19  is  denoted  by  line  A-B. 

Fig.  17c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  17d.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  17e.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  17f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  17g.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  17h.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  17i.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  123013,  at 
8.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  18.  Vertical  cross-section  of  reflectivity  and  wind.  View  is  toward  the  WNW  along 
the  cross-section  (A-B)  indicated  in  Fig.  17b. 

Fig.  19a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 

Fig.  19b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line.  Line  A-B  depicts  the  position 
of  the  cross-sectional  plot  depicted  in  Fig.  21. 

Fig.  19c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radax  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  19d.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radax  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  19e.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  19f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 


Fig.  19g.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  19h.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124913,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  20.  Vertical  cross-section  of  reflectivity  and  wind.  View  is  toward  the  WNW  along 
the  cross-section  (A-B)  indicated  in  Fig.  19b. 

Fig.  21a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  132902,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 

Fig.  21b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  132902,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line.  Line  A-B  denotes  the  position 
of  the  cross-sectional  plot  in  Fig.  22. 

Fig.  21c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  132902,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 

Fig.  22.  Vertical  cross-section  of  reflectivity  and  wind.  View  is  toward  the  WNW  along 
the  cross-section  (A-B)  indicated  in  Fig.  21b. 

Fig.  23.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  140905,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 

Fig.  24.  RHI  scan  of  reflectivity  and  radial  velocity  from  the  NCAR  CP-4,  (a)  is  re¬ 
flectivity  and  (b)  is  radial  velocity,  both  at  10  degrees  azimuth  at  110039  LST. 
Scales  are  given  at  the  bottom  of  each  figure. 
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Fig.  25.  RHI  scan  of  reflectivity  and  radial  velocity  from  the  NCAR  CP-4,  (a)  is  reflec¬ 
tivity  and  (b)  radial  velocity  for  30  degrees  azimuth  and  110108  LST.  Range 
marks  on  the  velocity  display  are  in  increments  of  10  km.  Scales  are  given  at 
the  bottom  of  each  figure. 

Fig.  26.  RHI  scans  of  reflectivity  and  radial  velocity.  Radar  is  NCAR  CP-4,  (a)  is 
reflectivity  and  (b)  is  radial  velocity,  both  at  210  degrees  azimuth  at  125830 
LST.  Scales  are  given  at  the  bottom  of  each  figure.  Range  marks  on  the  velocity 
display  are  in  increments  of  10  km. 

Fig.  27a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4  at 
1052.  Elevation  angle  is  1.8°.  Scales  are  given  at  the  bottom  of  each  figure. 

Fig.  27b.  VAD  at  1052.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
axe  given  by  the  thin  line. 

Fig.  27c.  VAD  at  1052.  Range  is  10.1  km.  Raw  velocity  values  a.e  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  27d.  VAD  at  1052.  Range  is  19.8  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  28a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4  at 
1152.  Elevation  angle  is  1.8°.  Scales  are  given  at  the  bottom  of  each  figure. 

Fig.  28b.  VAD  at  1152.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  29a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4  at 
1230.  Elevation  angle  is  1.8°.  Scales  are  given  at  the  bottom  of  each  figure. 

Fig.  29b.  VAD  at  1230.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  29c.  VAD  at  1230.  Range  is  10.1  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  29d.  VAD  at  1230.  Range  is  19.8  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  30a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4  at 
1249.  Elevation  angle  is  1.8°.  Scales  are  given  at  the  bottom  of  each  figure. 

Fig.  30b.  VAD  at  1249.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  30c.  VAD  at  1249.  Range  is  10.1  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  30d.  VAD  at  1249.  Range  is  19.8  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  31a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4  at 
1329.  Elevation  angle  is  1.8°.  Scales  are  given  at  the  bottom  of  each  figure. 

Fig.  31b.  VAD  at  1329.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  31c.  VAD  at  1329.  Range  is  10.1  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  31d.  VAD  at  1329.  Range  is  19.8  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  31e.  VAD  at  1329.  Range  is  29.9  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 

Fig.  Al.  Terminal  velocities  from  the  adjusted  Joss  and  Waldvogel  (1970)  equations. 

Plotted  lines  are  done  for:  (A)  10  dBZ,  (B)  20  dBZ,  (C)  30  dBz,  (D)  40  dBZ, 
and  (E)  50  dBZ.  Negative  velocities  denote  downward  speed.  All  plots  have  a 
value  of  -2  ms-1  above  the  10  km  height  level. 
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Radar  Analysis  of  A  TAMEX  Frontal  System 
Chapter  1 
Introduction 


1.1  The  TAMEX  Experiment 


J  The  Taiwan  Area  Mesoscale  Experiment  (TAMEX)  was  conducted  from  1  May 
through  29  June  1987.  It  was  a  joint  United  States-  Republic  of  China  research 
program.  Its  purpose  was  to  gather  data  to  enhance  the  understanding  of  the 
mesoscale  dynamical  and  microphysical  processes  responsible  for  the  rains  that  often 
lead  to  flash  flooding  in  Taiwan.  During  the  13  intensive  observational  periods 
(IOP’s),  data  was  gathered  on  1)  the  circulations  associated  with  the  MEI-YU 
front,  which  is  a  climatologically  favored  region  of  low-level  convergence  located 
over  Taiwan  and  southern  China;  2)  the  mesoscale  convective  systems  that  develop 
within  the  front;  and  3)  the  effects  of  the  topography  of  Taiwan  on  the  mean  flow, 
the  MEI-YU  front  and  mesoscale  convective  systems.  Njy#  , 

The  MEI-YU  front  usually  forms  between  the  tropical  and  polar  air  masses 
in  the  late  spring  and  early  summer.  It  moves  slowly  southeastward  to  a  quasi¬ 
stationary  position  extending  from  southern  Japan  to  southern  China.  Although 
the  MEI-YU  front  is  characterized  by  a  smaller  horizontal  temperature  gradient 
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them  is  found  in  typical  polar  fronts,  it  is  accompanied  by  a  strong  moisture  gradi¬ 
ent  and  horizontal  wind  shear.  Synoptically,  a  southwesterly  low  level  jet  is  often 
observed  to  the  south  of  the  MEI-YU  front.  Convection  usually  accompanies  the 
front.  Organized  mesoscale  convective  elements  are  embedded  within  the  frontal 
circulation  and  convection,  predominately  on  the  warm  side  of  the  front.  The 
convective  systems  generally  move  east  to  west,  sometimes  interacting  with  island 
terrain  to  produce  heavy  rainfall  (Chen  and  Yu,  1988). 

1.1.1  Geography 

The  dominant  topographic  feature  is  the  Central  Mountain  Range,  which  is 
300  km  in  length,  60-90  km  wide,  and  has  an  average  height  of  3000  m.  To  the 
west  of  the  mountains  is  a  coastal  plain  that  extends  inland  25-50  km.  On  the  east 
side  of  the  island  the  coastal  plane  extends  inland  only  about  10  km. 

The  three  Doppler  radars  (C-band)  were  located  on  the  northwest  coast  of 
Taiwan.  The  Chinese  Civil  Aeronautical  Administration  (CCA A)  radar  was  located 
at  Chaing  Kai-Shek  International  Airport.  The  National  Center  for  Atmospheric 
Research  (NCAR)  CP-4  radar  was  located  approximately  68  km  south  of  the  CCAA 
radar,  at  a  site  known  as  Waipu  Beach.  The  third  Doppler  radar,  the  National 
Oceanic  Atmospheric  Administration  (NOAA)  TOGA  radar,  was  located  at  the 
CCK  Air  Force  Base,  just  west  of  Taichung  and  approximately  45  km  south  of 
the  NCAR  radar.  Placed  in  these  areas,  this  radar  network  permitted  good  dual 
(sometimes  triple)  Doppler  coverage  of  many  of  the  IOP  phenomena  (Fig.  1). 


Fig.  1.  Doppler  radar  network  for  TAM  EX.  Cross  hatched  areas  denote  elevations 
greater  than  250  m,  and  stippled  areas  denote  elevations  above  1000  m.  From 
Cunning  (1988). 
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1.1.2  Synoptic  Setting 

The  focus  of  this  study  is  on  IOP  # 8 ,  which  occurred  from  1400  LST  7  June 
to  0800  LST,  9  June  (0000  GMT  corresponds  to  0800  LST).  The  observations  from 
this  period  were  selected  for  study  because  of  the  apparent  good  quality  of  Doppler 
radar  coverage  and  location  of  the  storms  with  respect  to  the  island.  The  radar 
data  analyzed  was  recorded  on  8  June  from  0800  LST  to  1630  LST.  The  objective 
of  this  IOP  was  to  observe  the  structure  of  the  convection  developing  over  the 
Central  Mountain  Range  (CMR)  of  Taiwan,  and  a  frontal  passage  over  the  island. 
It  was  this  frontal  area  that  was  studied.  The  ground-based  Doppler  radars  first 
observed  intense  convection  developing  over  the  mountains  on  7  June.  That  night, 
a  prefrontal  band  of  convection  developed  over  the  straits  and  moved  over  northern 
Taiwan.  During  the  second  day  (8  June),  the  Doppler  radars  collected  data  on  a 
second  band  of  convection  as  it  passed  over  the  island.  This  data  from  the  second 
band  is  the  data  that  was  analyzed  to  explore  the  similarities  and/or  differences 
with  mid-latitude  convection. 

A  surface  cold  front  on  7  June  was  over  mainland  China  moving  eastward 
towards  Taiwan,  with  a  weak  cloud  band  behind  the  front  and  scattered  embedded 
convection  within  the  band.  The  trough  was  expected  to  move  through  the  TAMEX 
area  within  the  next  36-48  hours.  There  was  some  weak  convection,  oriented  in  a 
north-south  band,  off  the  east  coast  of  Taiwan.  This  convection  band  may  have 
resulted  from  convergence  forcing,  created  by  the  outflow  from  earlier  convection 
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over  the  mountains  to  the  west.  The  local  soundings  reflected  a  slight  increase  in 
stability,  despite  the  fact  that  they  were  also  becoming  more  moist.  The  forecast 
was  for  the  front  to  move  into  the  Taiwan  area  at  about  midnight,  and  for  convection 
to  form  along  the  front  (Figs.  2  and  3). 

About  0100  LST,  8  June,  showers  developed  in  the  straits.  These  became 
better  organized,  moving  east-northeast,  and  by  0300  LST  the  line  was  in  the 
dual-Doppler  radar  area.  By  0500  LST,  the  line  moved  onshore  over  northern 
Taiwan,  although  most  of  the  line  missed  the  island.  A  0G00  LST  mesoscale  anal¬ 
ysis  by  Trier  ei  al.  (1989)  (See  Fig.  4)  revealed  light  winds  over  most  of  the  island 
with  weak  offshore  flow  observed  at  many  of  the  reporting  stations.  This  was  a 
typtVal  nocturnal  flow  pattern  observed  during  TAMEX,  during  undisturbed  condi- 
cions.  During  the  daylight  hours,  the  flow  tended  to  reverse  to  an  onshore  direction. 
For  undisturbed  conditions,  a  land-sea  breeze  circulation  is  typical.  At  this  time 
(0600  LST)  another  area  of  convection  moved  off  the  Chinese  mainland.  This  sec¬ 
ond  complex  of  convection  was  located  close  to  the  surface  front.  Prior  to  the  front 
reaching  Taiwan,  a  surface  pressure  ridge  was  observed  over  the  southwestern  side 
of  Taiwan,  with  a  lower  pressure  observed  on  the  eastern  side.  This  pattern  of  a 
windward  ridge  and  lee  side  trough  is  a  typical  pattern  according  to  Wang  (1986). 
This  pressure  pattern  sets  up  a  weak  southwesterly  flow  from  surface  to  4  km  over 
the  island.  As  shown  in  Fig.  5,  the  500  mb  trough  was  still  located  west  of  Taiwan. 
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Fig.  2.  Synoptic  scale  surface  analysis  for  TAMEX  on  1  June.  From  Cunning  (1988). 


Fig.  3.  Synoptic  scale  upper  air  analysis  for  TAMEX  on  7  June.  From  Cunning 
(1988). 
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Fig.  4-  Sea-level  pressure  with  1-mb  contour  intervals.  Station  model  symbols  have 
conventional  meanings  at  0600  LST  for  8  June.  From  Trier  et  al.  (1989). 
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The  forecast  called  for  the  front  to  continue  to  move  eastward  over  the  island. 
At  1200  LST,  the  surface  ridge  and  trough  pattern  persisted  over  much  of  the  island, 
except  over  the  northern  part,  where  analysis  at  this  time  showed  the  front  in  the 
Taiwan  straits  and  over  the  northwest  portion  of  the  island  (Fig.  6).  A  long  line 
of  frontal  convection  was  observed  with  maximum  reflectivity  values  of  40-45  dBZ. 
Trier  et  al.  (1989)  observed  that  this  front  was  characterized  by  sharp  temperature 
and  moisfure  contrasts.  Winds  observed  behind  the  front  shifted  from  southerly  to 
northeasterly,  paralleling  the  coast.  On  the  northeast  side  the  flow  turned  onshore. 
The  frontal  convection  continued  to  move  onshore  over  the  next  few  hours.  The 
convection  was  quite  extensive,  as  it  covered  the  area  from  north  of  the  island  down 
to  central  and  southern  Taiwan.  During  this  time  the  winds  in  the  higher  terrain 
of  the  Central  Mountain  Range  remained  predominately  southwesterly,  indicating 
that  the  front  was  a  shallow  phenomenon.  Lower  pressure  replaced  the  pressure 
ridge  on  the  windward  side.  The  flow  pattern  suggests  that  the  low-level  northerly 
flow  split,  and  was  being  channeled  around  the  Central  Mountain  Range  (Trier  et 
al.,  1989).  By  1800  LST,  most  of  the  convection  ended  in  the  dual-Doppler  area, 
although  some  convection  still  moved  onshore  over  southern  Taiwan,  as  shown  in 
Fig.  7. 

Stations  on  the  western  side  of  the  island  reported  that  after  the  front  had 
passed,  dewpoints  dropped  7°C,  and  pressure  rose  1.7  mb,  signaling  the  onset  of 
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convenUonat  meamng,  at  1100  LST  for  8  June.  From  Trier  et  al  (1989) 
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Fig.  7.  Sea-level  pressure  with  1-mb  contour  intervals.  Station  model  symbols  have 
conventional  meanings  at  1800  LST  for  8  June.  From  Trier  et  al.  (1989). 
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the  cold  air  mass.  Pressure  values  remained  stable  through  the  afternoon  and  then 
rose  significantly  in  the  evening  hours.  This  rise  may  have  been  the  result  of  an 
increase  in  the  depth  of  the  cold  air  accompanying  the  northerly  flow. 

The  eastern  seaboard  of  the  island  experienced  a  different  reaction  to  the  frontal 
passage.  In  the  east,  surface  heating  created  onshore  flow  by  late  morning,  causing 
the  flow  to  move  back  to  the  northeast  by  early  afternoon.  In  the  evening  hours,  the 
wind  pattern  returned  to  a  weak  offshore  flow  (even  after  frontal  passage).  This  was 
in  contrast  to  flow  on  the  west  coast  where  strong  northeasterly  flow  was  observed. 
The  dewpoints  on  the  east  coast  dropped  only  1°C.  The  east  coast  surface  pressure 
increase  was  the  only  real  indication  that  the  weather  had  changed.  Thus,  the 
strong  diurnal  forcing  masked  the  east  coast  frontal  passage.  Trier  et  al.  (1989) 
observed  that  the  wind  speed  behind  the  front  was  about  12  ms"1 .  At  0000  LST, 
9  June,  the  Doppler  radars  ceased  operations.  By  0800  LST,  9  June,  the  surface 
cold  front  moved  east  of  Taiwan  and  the  island  became  dominated  by  a  large  high 
pressure  system  (Fig.  8)  and  IOP  #8  ended. 

1.2  Tropical  Mesoscaie  Convective  Life  Cycle 

Rotunno  et  al.  (1988)  used  model  simulations  to  deduce  the  structure  and 
evolution  of  squall  lines.  In  accord  with  observations,  they  found  that  squall  lines 
are  fed  by  air  of  high  equivalent  potential  temperature,  8e ,  and  that  air  at  the  rear 


Fig.  8.  Synoptic  scale  surface  analysis  for  TAMEX  on  9  June.  From  Cunning  (1988). 
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of  the  line  is  cooler  and  has  a  lower  9e  than  the  air  in  front  of  it.  They  used  two- 
and  three-dimensional  models  to  examine  whether  a  squall  line  is  a  system  of  long- 
lived  cells,  or  a  long-lived  system  composed  of  ordinary  cells.  Their  experiments 
revealed  that  both  types  of  squall  lines  can  form.  The  type  that  forms  depends 
on  the  magnitude  and  orientation  of  the  environmental  shear  to  the  squall  line. 
Shallow  shear,  oriented  perpendicular  to  the  squall  line,  develops  a  long-lived  line 
with  short-lived  cells,  suggesting  an  interaction  between  cells.  Cold  surface  outflow 
and  low-level  shear  creates  much  deeper  lifting  than  is  possible  without  the  shear, 
allowing  for  new  cell  formation  as  the  old  cells  decay.  With  deep  shear,  a  line  of 
long-lived  cells  form.  The  environment  of  the  slow  mover  is  characterized  by  little 
or  no  shear,  and  a  significantly  more  moist  layer  above  the  boundary  layer  than 
that  of  the  fast  mover.  Rotunno  et  al.  (1988)  showed  that  the  value  of  the  shear 
necessary  to  maintain  a  long-lived  state  depends  on  the  pr  operties  of  the  cold  pool. 
The  cold  pool,  in  turn,  is  dependent  on  the  thermodynamic  sounding  and  the  type 
of  the  convection  created.  As  the  temperature  in  the  cold  pool  decreases,  it  begins 
to  overcome  the  low-level  shear.  It  then  starts  to  surge  outward  underneath  the 
updraft  in  the  downshear  direction.  This  creates  an  upshear  tilt;  and  it  is  this 
upshear  tilt  that  signals  the  dominance  of  the  cold  pool  circulation  which  leads  to 
less  intense  convection. 

Tropical  convection  is  often  classified  as  a  Mesoscale  Convective  System  (MCS) 
because  tropical  convection  shares  a  similar  structure  and  a  common  dynamical  ba- 
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sis  with  mid-latitude  convection.  Leary  and  Houze  (1979)  studied  tropical  convec¬ 
tion  using  GARP  Atlantic  Tropical  Experiment  (GATE)  data.  In  this  study  they 
examined  six  mesoscale  cloud  clusters  that  occurred  on  5  September  1974.  The 
focus  of  the  study  was  on  the  structure  and  evolution  of  the  convection.  The  con¬ 
vective  cycle  is  sub-divided  into  four  stages:  the  formative  stage,  the  intensifying 
stage,  the  mature  stage,  and  the  dissipating  stage. 

1.2.1  The  Formative  Stage 

In  the  formative  stage  the  tropical  Mesoscale  Convective  System  (MCS)  con¬ 
sists  of  a  line  of  isolated  echoes.  When  a  MCS  organizes  into  a  linear  system  it  is 
called  a  squall  line.  There  are  special  conditions  that  are  believed  to  be  required 
for  the  formation  of  squall  lines.  Most  important  is  an  unstable  air  mass,  that  is, 
an  air  mass  that  is  warm  and  moist  at  low  levels  and  relatively  colder  aloft.  A 
dryline  above  the  warm  moist  layer  enhances  the  vigor  of  the  squall  line.  Precip¬ 
itation  that  evaporates  in  this  dry  region  cools  the  air  and  increases  the  negative 
buoyancy,  which  in  turn  strengthens  the  downdraft  and  outflow.  Once  the  squall 
line  has  formed,  it  will  be  able  to  maintain  itself  through  convergence  along  the 
outflow  boundary.  This  causes  strong  ascent  into  existing  updraft  areas,  or  initiates 
convection  ahead  of  the  storms  that  already  exist  (Hane,  1986). 

Most  convection  forms  in  a  conditionally  unstable  environment  with  a  trig¬ 
gering  mechanism.  An  unstable  environment  is  insufficient  to  initiate  squall  line 
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convection;  a  triggering  mechanism  is  required.  An  example  of  a  triggering  mech¬ 
anism  is  a  cold  front.  The  triggering  mechanism  can  be  a  cold  front  or  a  dryline. 
The  intersection  of  a  front  and  a  dryline  can  produce  frontogenesis,  as  examined 
by  Kessinger  and  Bluestein  (1979)  and  Bluestein  and  Berry  (1982).  Newton  and 
Fankhauser  (1964)  examined  a  method  in  which  a  squall  line  grew  by  adding  new 
cells  to  the  southern  edge  of  the  line.  This  building  occurs  either  as  a  result  of  an 
upper  level  disturbance  moving  south  over  a  surface  discontinuity,  or  because  the 
outflow  from  a  storm  interacts  with  inflow  at  the  southern  end  of  the  line.  This 
forms  a  convergence  zone,  and  thus  an  area  conducive  to  new  eonvectional  growth. 

Leary  and  Houze  (1979)  stated  that  echoes  viewed  in  a  vertical  cross-section 
have  vertically-oriented  contours  of  reflectivity.  This  initial  formation  is  typical 
of  situations  that  develop  into  areas  of  high  intensity,  as  well  as  lines  of  echoes 
that  dissipate  without  any  further  development.  Thus,  the  factors  that  control  the 
development  are  not  structural,  but  dynamical.  The  lines  that  dissipate  without  an 
increase  in  development  are  found  to  be  aligned  nearly  parallel  to,  and  move  in  the 
same  direction  as,  the  surface  wind.  The  lines  of  convection  that  are  observed  to 
intensify  are  found  to  be  aligned  more  nearly  perpendicular  with,  and  move  in  the 
opposite  direction  of,  the  surface  wind.  The  low-level  outflow  from  the  convective 
scale  downdrafts  in  these  echoes  are  aligned  perpendicular  to  the  surface  wind.  This 
causes  a  modification  of  the  local  low-level  convergence  pattern,  making  it  conducive 
for  new  cell  formation  in  the  upwind  direction.  Further,  if  the  convective  downdraft 


is  assumed  to  propagate  outward  in  a  radial  fashion,  then  low  level  convergence 
(and  thus,  new  cell  formation)  is  favored  between  the  individual  echoes  in  the  line. 

1.2.2  The  Intensifying  Stage 

In  the  intensifying  stage,  development  in  the  upwind  direction  contributes  to 
the  propagation  speed  of  the  mesoscale  precipitation  region.  New  cells  form  along 
the  leading  edge  of  precipitation  before  the  older  cells  at  the  rear  can  dissipate. 
This  broadens  the  mesoscale  line.  When  new  cells  are  generated  between  the  older 
cells,  they  merge  with  each  other.  The  resulting  formation  appears  as  a  large  echo, 
with  many  intense  convective  cores.  These  echoes  grow  in  height,  and  the  upper 
level  winds  carry  the  cloud  and  precipitation  particles  downwind,  forming  the  cirrus 
shield,  or  anvil. 

1.2.3  The  Mature  Stage 

Leary  and  Houze  (1979)  state  the  mature  stage  resembles  the  intensifying  stage 
in  many  respects.  In  the  mature  stage,  the  echoes  exhibit  extensive  development 
both  in  surface  precipitation  pattern  and  in  the  vertical  cross  section.  The  deep 
convective  cells  at  the  leading  edge  overlie  the  most  intense  surface  precipitation. 
Upper  level  winds  continue  to  carry  cloud  and  precipitation  particles  downwind. 
The  major  difference  in  the  structure  between  the  mature  aim  the  intensifying  stage 
is  that  in  the  mature  stage  an  extensive  region  of  horizontally  uniform  precipitation 
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is  to  the  rear  of  the  most  convective  cells.  This  region  to  the  rear  is  composed  of 
cells  from  the  leading  edge  of  the  line  that  have  weakened  as  they  blend  together, 
and  “moved”  to  the  rear  of  the  precipitation  region.  In  this  horizontally  uniform 
region,  often  a  bright  band  of  radar  reflectivity  can  often  be  seen.  This  bright 
band  is  located  in  the  melting  layer,  just  below  the  0°C  isotherm.  This  horizontally 
stratified  reflectivity  pattern  is  inconsistent  with  the  patterns  found  in  areas  of 
convective  scale  up  and  downdrafts  (Leary  and  Houze,  1979). 

1.2.4  The  Dissipating  Stage 

The  transition  from  the  mature  stage  to  the  dissipating  stage  is  a  gradual  one. 
Intense  convective  cells  continue  to  form  along  the  leading  edge,  but  they  become 
progressively  shallower,  and  eventually  cease  to  form  at  all.  Then  the  large  area  of 
uniform  precipitation  begins  to  slowly  dissipate  until  only  a  few  fragments  of  the 
original  mass  remains  (Leary  and  Houze,  1979). 

1.3  Echo  Scales  and  Sizes 

Houze  and  Cheng  (1977)  analyzed  the  echoes  observed  during  the  GATE  ex¬ 
periment.  Although  their  study  did  not  examine  frontal  systems  specifically,  it  did 
describe  the  radar  echoes  in  terms  of  area,  heights,  duration,  intensity,  motion,  and 
the  orientation  of  the  echo  lines.  Their  study  was  based  on  the  analysis  of  radar 
data,  but  it  did  not  involve  radar  intercomparisons,  calibration  studies  or  extensive 
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data  processing.  It  was  intended  to  serve  as  a  predecessor  for  more  quantitative 
analysis. 

The  horizontal  area  of  precipitation  was  divided  into  a  series  of  size  scales.  The 
scales  were  denoted:  the  B/C  scale  for  103-105  km2,  the  C  scale  for  102-103  km2, 
and  the  D  scale  for  areas  1-102  km2 ,  or  less.  The  population  of  radar  echoes  recorded 
in  GATE  were  predominately  in  the  D  scale  range  (67%)  and  the  C  scale  accounted 
for  28%.  These  echoes,  however,  covered  only  21%  of  the  total  area  of  echoes.  Thus, 
only  a  relatively  few  large  precipitating  areas  are  found  in  tropical  cloud  populations, 
but  these  take  on  importance  when  their  areal  extent  is  considered.  The  average 
height  for  scale  echoes  D  was  4.3  km;  for  C  scale  echoes,  the  average  height  was  5.8 
km,  and  for  B/C  scale  echoes,  the  average  height  was  9.9  km.  The  average  height 
of  the  GATE  echoes  was  4.8  kin.  Echo  duration  was  positively  correlated  with  the 
echo  area  coverage.  The  overall  mean  echo  duration  was  88  minutes.  Most  of  the 
echoes  lasted  less  than  50  minutes.  Echoes  with  high  intensity  cores  were  most 
often  observed  in  the  B/C  regime.  Echoes  that  exceeded  104  km2  in  area  contained 
as  many  high  intensity  cores  as  moderate  and  low-level  intensity  cores.  Higher 
intensity  cores  were  more  often  found  in  the  B/C  scale,  indicating  that  these  large 
echoes  form  in  an  environment  favorable  for  the  development  of  intense  cores,  or 
that  the  more  intense  cores  created  surrounding  areas  of  widespread  precipitation. 
Houze  and  Cheng  (1977)  calculated  that  there  was  an  average  of  one  echo  core  for 
every  187  km2  of  echo  area. 
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James  and  Browning  (1979)  used  Doppler  radar  to  study  15  examples  of  frontal 
convection.  The  Doppler  radar  data  used  in  this  study  was  collected  by  Meteorolog¬ 
ical  Office  Radar  Research  Laboratory  during  the  years  1968-1977.  The  convection 
was  often  broken  into  cells  with  the  average  length  of  an  element  22  km  and  the 
average  gap  between  cells  9  km.  Occasionally,  the  features  ranged  in  size  from  a 
few  kilometers  to  100  km.  The  gaps  between  cells  also  varied  in  length.  They 
were  usually  a  few  kilometers  to  tens  of  kilometers  wide,  but  in  one  case  they  were 
observed  to  be  95  km  wide.  The  lifetime  of  the  elements  was  also  found  to  vary, 
ranging  from  10  minutes  to  4  hours.  However,  James  and  Browning  (1979)  found 
that  linear  convection  can  be  up  to  200  km  long,  with  no  elemental  structures  pre¬ 
sent.  In  this  situation,  areas  of  higher  echo  intensity  can  be  observed  within  the 
line.  They  found  no  distinguishing  thermodynamic  or  dynamic  characteristics  that 
would  account  for  the  difference  in  the  elemental  and  non-elemental  structures.  The 
convection  in  the  cases  with  no  elements  was  slower  moving  and  exhibited  a  slightly 
greater  vector  wind  and  temperature  change  across  the  frontal  zone. 


1.4  Echo  Morphology 


Continuity  and  geometric  relationship  of  the  echoes,  and  the  persistence  of 
precipitation  are  some  of  the  characteristics  of  a  squall  line.  Squall  lines  tend  to 
favor  formation  in  areas  of  quasi-linear  forcing  in  the  warm  sector. 
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Kessler  and  Wexler  (1959)  used  an  incoherent  radar  to  observe  a  passing  cold 
front  in  New  England  in  one  of  the  first  observations  of  a  cold  front  with  a  radar. 
The  front  was  characterized  by  strong  gusty  winds  and  heavy  rain  showers  with  a 
squaili  line  ahead  of  the  front.  Convection  in  the  squall  line  was  limited  to  a  height 
of  4  km,  whereas  the  convection  closer  to  the  front  reached  heights  up  to  10  km. 
The  most  intense  convection  was  below  4  km,  where  temperatures  were  greater  than 
0°C.  The  cold  front  had  strong  vertical  velocities  at  low  levels,  and  unusually  heavy 
precipitation. 

Browning  and  Harrold  (1970)  used  a  Doppler  radar  to  study  two  cold  fronts 
that  passed  through  Great  Britain.  They  studied  the  two-dimensionality  of  the 
frontal  structure  over  a  25  km  length.  Intense  echoes,  extending  to  3  km  in  the 
vertical,  were  observed  at  the  leading  edge.  This  is  comparable  to  the  heights 
found  by  Kessler  and  Wexler  (1959).  Behind  the  front,  the  tops  of  the  echoes 
extended  up  to  a  height  of  6  km.  Behind  the  line  of  most  active  convection,  a 
500  m  thick  layer  of  warm  air  ascended  at  a  small  angle  relative  to  the  surface. 
The  base  of  this  warm-air  layer  was  1.5  km  height  immediately  behind  the  front, 
ascending  to  a  height  of  4  km,  100  km  behind  the  front.  Radar  data  showed  that  the 
line  was  devoid  of  pronounced  cellular  structure.  Cellular  patterns  are  commonly 
associated  with  cumuloform  convection.  Thus,  Browning  and  Harrold  (1970)  stated 
that,  to  a  first  approximation,  the  flow  in  the  surface  cold  front  could  be  treated 
as  two-dimensional.  This  assumption  was  further  justified  when  Browning  and 
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Harrold  examined  8  vertical  sections  that  were  taken  over  a  25  minute  period  along 
a  24  km  section  of  the  front.  Each  of  these  vertical  sections  had  similar  airflow 
patterns,  supporting  the  two-dimensional  approximation.  An  interesting  feature  of 
the  observed  echo  pattern  was  the  occasional  appearance  of  a  vault.  A  vault  is  a 
region  of  relatively  (compared  to  its  surroundings)  weak  echo  within  the  updraft 
core.  A  vault,  up  to  this  point,  had  only  been  observed  in  severe  local  storms. 
In  this  case,  the  vault  was  1-2  km  wide  and  the  echo  was  15  dB  less  than  the 
surrounding  echo  area. 

Hobbs  et  al.  (1980)  examined  a  passing  fronted  system  that  was  125  km  wide 
and  composed  of  five  mesoscale  rainbands.  The  data  was  recorded  by  NCAR’s 
CP-3  radar  on  17  November  1976,  during  the  University  of  Washington’s  Cyclonic 
Extratropical  Storms  (CYCLES)  Project.  Four  of  these  rainbands  were  several  tens 
of  kilometers  wide,  and  the  other  band,  which  was  coincidental  with  the  surface 
front,  was  only  4  km  wide.  These  bands  contained  precipitation  cores  that  varied  in 
size  from  25  to  300  km2 .  The  wide  cold-frontal  rainbands  had  irregular  shapes  and 
were  randomly  positioned.  The  cores  associated  with  the  narrow  cold-frontal  bands 
were  elliptical.  The  height  of  the  core  tops  were  determined  by  the  Doppler  radar 
observations.  The  tops  of  the  wide  rainbands  extended  to  5-6  km.  The  tops  of  the 
narrow  band  were  between  1.5-4. 5  km  in  height.  Doppler  measurements  were  used 
to  deduce  the  wind  directions  throughout  the  frontal  system.  Near  the  leading  edge 
of  the  frontal  system,  a  one  km  thick  overhanging  cloud  of  precipitation  extended 
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approximately  30  km  in  front  of  the  system.  Another  overhanging  area  5  km  thick 
extended  90  km  behind  the  main  body  of  the  frontal  system.  Precipitation  from 
this  rearward  extending  overhang  evaporated  prior  to  reaching  the  ground. 

In  general,  the  observations  have  shown  that  the  echoes  exhibit  certain  charac¬ 
teristics:  a  line  of  intense  echoes,  about  4  km  high  at  the  leading  edge  of  the  echo 
region,  followed  by  an  area  of  echoes  that  increase  in  height,  corresponding  to  the 
front  aloft.  Finally,  the  examination  of  cross-sections  revealed  a  two-dimensional 
character  of  the  system. 

1.5  Flow  Patterns,  Convergent  and  Divergent  Patterns 

The  understanding  of  dynamical  forces  and  interactions  has  been  greatly  ad¬ 
vanced  by  three-dimensional  time-dependent  numerical  models.  Coupled  with  this 
have  been  the  advances  in  dual-  and  triple-Doppler  radar  analysis  to  illuminate  the 
internal  storm  structure. 

Bergeron  (1954)  and  Pedgley  (1962)  studied  the  flow  patterns  within  a  squall 
line.  They  observed  that  the  region  of  stratiform  rain  ascended  at  high  levels  and 
descended  at  low  levels.  An  area  of  strong  convection  was  found  at  the  leading 
edge  of  the  system.  Smull  and  Houze  (1985,  1987)  used  Doppler  radar  data  to 
show  that  updraft  flow  exists  ahead  of  the  line’s  outflow,  and  that  this  updraft  tilts 
backwards  into  the  line.  Another  characteristic  that  they  found  was  a  region  of 
back-to-front  flow  near  the  post-line  rainbund  at  3  8  lun  height  level.  This  rear 
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inflow  is  not  commonly  seen  with  tropical  squall  lines.  This  is  most  likely  because 
the  tropics  lack  the  strong  instability,  strong  cold  pools,  and  strong  shear  needed 
for  maintenance  of  this  flow  (Ray,  1988). 

Kessler  and  Wexler  (1959)  observed  that  by  increasing  the  gain,  they  were  able 
to  observe  gaps  of  a  few  kilometers  between  the  echoes  from  the  heavy  rain  and  the 
lighter  precipitation  on  the  upper  side  of  the  front.  This  suggested  that  the  updraft 
currents  were  associated  with  adjacent  compensating  downdrafts. 

Browning  and  Harrold  (1970)  analyzed  a  cold  front  by  tracing  the  pattern 
of  motion  using  the  wet  bulb  potential  temperature  as  a  tracer  of  airflow.  They 
observed  a  narrow  tongue  of  warm  air  just  ahead  of  the  front.  This  warm  air  was 
advected  (relative  to  system  movement)  by  strong  geostrophic  flow,  moving  parallel 
to  the  front.  It  was  this  prefrontal  warm  tongue,  ascending  at  the  front,  that  was 
responsible  for  the  growth  of  frontal  precipitation.  They  treated  the  flow  in  this 
warm  tongue  in  two  parts.  The  first  part  was  for  the  1-3  km  height  flow,  which  was 
parallel  to  the  front.  Below  the  1  km  height,  winds  backed  due  to  surface  friction, 
and  a  component  was  observed  toward  the  cold  front.  The  overall  flow  pattern  in 
the  warm  tongue  showed  mass  flux  toward  the  front  in  the  lowest  kilometer  of  the 
frontal  zone.  This  was  fed  by  subsiding  air  from  the  flow  parallel  to  the  front.  The 
flux  toward  the  cold  front  was  seen  to  rise  when  it  entered  the  line  of  convection 
accompanying  the  surface  cold  front. 
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The  warm  air  flow  parallel  to  the  front  had  a  component  of  velocity  of  24  ms-1 
parallel  to  the  front,  and  a  component  of  2  ms-1  toward  it.  Below  1  km,  flow  into 
the  front  was  in  a  region  where  convergence  values  were  10-2s-1.  The  updraft  flow 
was  seen  tilted  60°  from  the  horizontal,  reaching  a  maximum  speed  of  8  ms-1  at 
1.2  km  in  height.  The  updraft  extended  to  3  km  height.  Below  this,  the  horizontal 
divergence  values  were  as  large  as  9  x  10-3s-1.  Some  of  this  divergent  air  flowed 
into  the  area  ahead  of  the  front,  but  most  of  it  departed  the  surface  frontal  updraft 
area  in  the  rear  of  the  front  between  the  heights  of  1  and  3  km.  This  air  then  began 
to  ascend  on  a  gentle  slope,  as  previously  discussed.  Downdraft  flows  were  almost 
always  less  than  1  ms-1  in  the  cold  air,  which  is  consistent  with  the  absence  of 
pronounced  surface  gusts,  or  humidity  drops,  both  of  which  are  commonly  observed 
in  squall  lines. 

Radiosonde  data  taken  prior  to  frontal  passage  revealed  a  neutral  lapse  rate 
at  almost  all  levels.  The  Doppler  observed  streamlines  showed  that  the  cold  air 
did  not  overrun  the  warm  air  mass.  Instead,  the  warm  air  rose  over  the  advancing 
cold  air.  Thus,  the  convection  line  was  produced  by  forcing,  rather  than  by  free 
convection.  This  forced  convection  viewpoint  is  supported  by  the  absence  of  cellular 
structure  at  the  surface  cold  front.  Almost  the  entire  updraft  at  the  surface  was  the 
result  of  frictional  convergence.  The  large  flux  of  air,  and  hence  the  precipitation, 
were  associated  with  the  large  values  of  flow  parallel  the  front  at  the  leading  edge 
(Browning  and  Harrold,  1970). 
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Behind  the  nearly  vertical  ascent  at  the  surface,  the  warm  air  gradually  as¬ 
cended  towards  the  rear  of  the  system.  This  warm  air  was  found  in  a  shallow  layer 
between  the  advancing  cold  air  wedg.;  and  the  stable  layer  at  3  km  height.  The 
circulation  in  this  area  was  well  marked,  with  updraft  values  in  the  warm  air  of 
20  cm  s-1.  The  descent  in  the  cold  air  was  also  evident,  but  it  was  weaker  com¬ 
pared  to  the  ascending  warm  air.  This  was  not  just  a  simple  case  of  warm  air  riding 
up  and  over  the  advancing  cold  wedge,  because  the  transition  area  from  ascent 
to  descent  was  within  the  warm  air  instead  of  the  frontal  zone'.  The  cloud  base, 
assumed  to  be  at  the  level  of  this  transition,  was  in  the  warm  air.  The  cold  air 
descended,  with  the  strongest  area  of  descent  in  the  frontal  zone  itself.  The  down¬ 
ward  velocities  of  the  air  in  the  frontal  zone  was  comparable  to  the  ascending  warm 
air  (20  cm  s-1).  The  maximum  wind  speed  toward  the  surface  front  was  8  ms-1. 
Browning  and  Harrold  (1970)  further  noted  that  although  the  cold  frontal  zone  was 
characterized  by  an  area  of  descent,  there  was  an  apparent  downward  extension  of 
the  region  of  ascent  into  the  frontal  zone.  This  feature  may  have  been  associated 
with  an  undulation  in  the  frontal  zone  boundaries  such  that  air  remained  within 
this  zone.  It  also  may  have  been  the  result  of  a  small  flux  of  air  upward  through 
the  frontal  zone.  Doppler  radar  resolution  of  this  region  showed  that  the  airflow 
pattern  was  fluctuating  rapidly  on  a  scale  of  kilometers,  indicating  that  this  was  an 
area  of  strong  mixing.  Mixing  is  consistent  with  earlier  studies.  In  the  lowest  few 
hundred  meters,  the  cold  air  was  frictionally  retarded  to  the  point  that  a  component 
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was  seen  in  the  direction  away  from  the  cold  front.  The  depth  of  this  rearward  flow 
increased  from  near  zero  at  the  front  to  1  km  at  80  km  behind  the  front. 

Parsons  et  al.  (1987)  deduced  the  pressure  fields  in  a  squall  line.  A  high  pres¬ 
sure  area  was  found  below  the  maximum  updraft,  with  pressure  increasing  at  low 
levels  when  going  from  the  warm  sector  to  the  cold  sector,  and  a  low  pressure  above 
the  cold  air  with  a  pressure  minimum  above  the  front.  It  was  this  vertically-oriented 
pressure  gradient  that  created  the  updrafts  in  the  absence  of  potential  instability. 
Examination  of  the  buoyancy  fields  showed  no  evidence  for  buoyancy-created  up¬ 
drafts.  Thus  the  theory  of  forced  convergence  and  vertical  pressure  gradient  forced 
updraft  was  further  supported  (Parsons  et  al.,  1987). 

Moore  and  Heggli  (1979)  studied  a  winter  squall  line  and  cold  front  that  passed 
through  California  on  5  February  1978.  A  low-level  jet  of  25  ms-1  was  located  just 
ahead  of  the  front.  Satellite  cloud  band  tracking  showed  that  the  front  moved  at 
about  30  ms-1  over  the  Pacific  ocean.  Intense  convection  was  observed  by  satellite 
2-3  hours  before  landfall.  A  general  area  of  subsidence  was  indicated  around  the 
squall  line.  Radar  analysis  of  the  front  showed  that  the  low-level  jet  was  converging 
with  the  southwesterly  winds  from  behind  the  front.  The  radar  recorded  reflectiv¬ 
ities  exceeding  50  dBZ.  The  size  of  the  higher  reflectivities  areas  extended  up  to 
2  km  wide.  The  reflectivities  lacked  the  characteristics  of  cellular  type  formation. 
Other  features  were  seen  on  the  1.1°  elevation  scan.  An  area  of  subsidence  was  seen 
proceeding  the  frontal  band.  This  area  extended  50  km  ahead  of  the  front.  The 
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frontal  position,  as  determined  by  the  velocity  discontinuity  line  on  the  Doppler 
radar,  was  placed  just  ahead  of  the  line  of  the  strongest  reflectivities. 

Soundings  taken  prior  to  frontal  passage  revealed  that  the  environment  was,  in 
general,  neutral  to  slightly  stable.  Doppler  measurements  showed  that  the  low-level 
jet  (1.2  km  in  height)  and  the  area  below  it  was  the  source  of  air  drawn  into  the 
front  as  it  passed.  Updraft  velocities  reached  17  ms-1  at  the  1.8  km  height.  One 
prefrontal  unstable  level  was  found  between  0.5  and  2.0  km  in  height.  By  comparing 
the  observed  velocities  with  the  velocities  that  could  have  been  created  by  positive 
buoyancy  (15  ms-1,  in  this  case),  Moore  and  Heggli  (1979)  concluded  that  the 
vertical  motions  were  due  to  boundary  layer  forcing,  and  not  hydrostatic  instability. 
Post-frontal  soundings  revealed  neutral  stability  below  1.7  km,  an  isothermal  layer 
(associated  with  the  upper  level  front)  from  1. 7-2.0  km,  and  stable  layers  above  2.0 
km.  Horizontal  divergence  was  seen  just  above  the  frontal  boundary.  It  was  also 
observed  that  as  the  front  moved  over  land  it  propagated  at  approximately  20  ms-1 , 
but  it  slowed  to  about  10  ms-1  as  it  encountered  the  mountains. 

Newton  (1950)  and  Newton  and  Newton  (1959)  showed  that  new  convection 
formed  in  the  convergence  zone  of  the  outflow  boundary.  They  stated  that  this  was 
maintained  by  the  transfer  of  westerly  momentum  downward.  Kessler  and  Wexler 
(1959)  also  observed  very  large  horizontal  convergences  confined  to  a  very  narrow 
zone  near  the  ground.  This  created  a  rapid  increase  in  the  vertical  velocities  with 
height  at  the  lower  levels.  The  radar  also  indicated  the  presence  of  strong  horizontal 
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divergence  above  the  frontal  zone  at  the  2.5  km  to  4.0  km  levels.  The  radar  de¬ 
tected  heavy  rain  that  maintained  its  position  along  the  front,  moving  at  the  same 
speed  as  the  front.  Kessler  and  Wexler  (1959)  also  observed  that  cell  development 
appeared  to  occur  on  the  forward  edge  of  the  heavy  rain  line,  and  that  dissipation 
occurred  on  the  back  edge.  The  convective  character  of  the  front’s  precipitation 
was  seen  to  strengthen  as  a  result  of  daytime  heating.  Strong  northwest  winds  were 
recorded  that  exceeded  the  winds  aloft,  as  measured  by  radiosonde.  Thus,  it  was 
determined  that  these  gusts  could  not  be  due  to  transport  of  momentum  through 
vertical  motions  alone.  Kessler  and  Wexler  theorized  that  the  gusts  were  probably 
due  to  localized  strong  pressure  gradients  near  the  front. 

Matejka  and  Houze  (1978)  analyzed  the  airflow  in  a  vertical  cross-section  of  the 
cold  front.  The  data  was  recorded  by  NCAR’s  CP-3  radar  on  17  November,  1976, 
during  the  University  of  Washington’s  Cyclonic  Extratropical  Storms  (CYCLES) 
Project.  To  investigate  the  airflow  within  the  cold  front,  they  applied  the  equation 
of  continuity  to  the  horizontal  component  of  the  velocity  in  a  single  vertical  plane, 
oriented  perpendicular  to  the  front.  The  vertical  velocities  were  calculated  assum¬ 
ing  that  no  divergence  occurred  along  the  front,  and  the  the  cross-front  divergence 
was  integrated  upward  in  one-third  km  intervals  from  the  surface.  This  was  done 
using  zero  for  vertical  velocity  at  the  surface.  The  analysis  of  the  cross-front  di¬ 
vergence  showed  that  the  strongest  convergence  (1.9  x  10~3s_1)  occurred  at  the 
surface  windshift  in  the  lowest  kilometer.  Weaker  convergence  straddled  the  cold 


31 


front  aloft  and  divergence  was  seen  above  and  ahead  of  the  front.  The  strongest 
area  of  divergence  was  at  a  height  of  1.5  km,  immediately  in  front  of  the  surface 
wind  shift.  The  pattern  of  vertical  velocity  revealed  two  distinct  regions  of  ascent. 
One  region,  intense  and  small,  was  found  in  the  area  of  convergence  associated  with 
the  wind  shift.  The  other  region  of  ascent  was  larger  and  was  located  above  the 
cold  frontal  surface  aloft.  The  airflow  in  the  warm  sector  moved  toward  the  front 
(in  a  relative  sense)  at  speeds  of  13  ms""1.  When  it  reached  the  front,  this  air  rose 
in  a  column  4  km  wide  at  1  ms-1.  When  the  air  reached  approximately  3  km  in 
height,  wind  carried  it  forward.  Thus,  it  is  possible  that  precipitation  from  this  area 
may  have  fallen  into  the  convection  at  the  leading  edge  of  the  system,  ahead  of  the 
surface  wind  shift.  A  0.5  ms-1  downdraft,  4-5  km  wide,  was  observed  just  ahead  of 
the  updraft.  The  base  of  this  downdraft  was  at  the  1.5  km  height,  which  coincided 
with  the  location  of  maximum  divergence.  The  updraft-downdraft  couplet  observed 
by  Matejka  and  Houze  is  typical  of  cumulonimbus  storms.  Thus,  they  were  able 
to  distinguish  smaller  scale  flows  within  the  broader  band  of  frontal  motions.  The 
broader  region  of  ascent  associated  with  the  cold  front  aloft  had  vertical  motions  of 
a  fraction  of  a  ms-1.  This  was  a  result  of  the  convergence  that  straddled  the  front 
aloft.  Matejka  and  Houze  (1978)  also  observed  that  this  ascent  was  not  uniform 
throughout  the  frontal  band,  and  in  areas  where  the  velocities  were  stronger  heav¬ 
ier  precipitation  was  observed.  Their  results  confirm  the  earlier  results  of  Browning 
and  Harrold  (1970)  showing  two  distinct  regions  of  ascent. 
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Carbone  and  Serafin  (1979)  studied  the  same  front  as  Moore  and  Heggli  (1979). 
Their  data  set  had  two  limitations-one  was  that  a  power  loss  prevented  data  acqui¬ 
sition  above  elevation  angles  of  15°.  The  other  limitation  was  that  surface  horizon 
prevented  two  of  the  radars  from  data  collection  below  0.9  km.  The  data  was  in¬ 
terpolated  to  a  Cartesian  grid  using  an  isotropic  Gaussian  weighting  function,  with 
a  300  m  half-weight  radius.  The  sphere  of  influence  extended  to  750  m  in  radius. 
The  grid  spacing  was  300  m  in  all  directions.  This  filter  radius  and  grid  spacing 
combination  resolves  motions  greater  than  1.2  km  in  wavelength.  Mass  continuity 
was  used  to  indirectly  obtain  vertical  motions.  An  upper  boundary  condition  of 
w  =  0  at  z  =  6.3  km  was  used.  This  value  was  in  line  with  the  hydrostatically  sta¬ 
ble  environment,  especially  when  one  takes  into  account  the  stratiform  reflectivity 
structure  and  the  weak  divergence  at  the  storm  top.  Downward  integration  of  the 
continuity  equation  was  the  final  step.  The  general  circulation  of  the  front  and  the 
squall  line  remained  steady,  until  tornadic  formation  began  to  occur. 

When  the  eddy  velocity  field  was  examined,  at  the  3.6  km  height,  strong  two- 
dimensional  divergence  was  observed.  Carbone  and  Serafin  (1979)  also  examined  the 
storm  relative  wind  field.  In  doing  this,  flow  parallel  to  the  main  band  was  seen.  This 
denoted  the  area  of  reflectivity  overhang.  Most  of  the  diverging  updraft  exits  to  the 
rear  of  the  line.  Two-dimensional  cross  sections  of  the  line  revealed  an  up/downdraft 
couplet  that  was  similar  to  those  found  in  multiple-Doppler  syntheses.  Carbone  and 
Serafin  (1979)  observed  that  the  vertical  updraft  was  4  km  wide,  and  reached  speeds 
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of  17  ms-1  at  1.8  km  height,  in  agreement  with  Moore  and  Heggli  (1979).  Further, 
Carbone  and  Serafin  (1979)  found  that  the  updraft  became  divergent  at  2.1  km, 
decreasing  in  speed,  and  reached  a  maximum  height  of  6.0  km.  The  majority  of  this 
air  exits  the  storm  at  the  rear  flank.  The  maximum  downdraft  speed  was  12  ms-1, 
and  was  found  immediately  behind  the  surface  convergence  line.  This  was  also  the 
area  of  the  highest  reflectivities  and  greatest  rainfall.  From  the  southwest,  below 
2  km,  westerly  momentum  air  was  seen  approaching  the  convergence  line.  Shallow 
post-frontal  clouds  were  seen  in  this  area.  The  average  vertical  velocity  was  3  ms-1. 

Hobbs  et  al.  (1980)  conducted  a  more  detailed  study  of  the  cold  front  examined 
by  Matejka  and  Houze  (1978).  Not  unexpectedly,  many  of  the  results  in  this  study 
are  in  good  agreement  with  those  from  the  earlier  work.  The  analysis  by  Hobbs 
et  al.  (1980)  examined  the  frontal  system  over  a  longer  time  period,  and,  thus,  it 
included  a  study  of  the  front  prior  to  landfall.  The  radar  was  used  to  obtain  the 
cross-band  component  of  the  horizontal  wind  in  vertical  planes  perpendicular  to 
the  rainbands.  Vertical  velocities  were  calculated  using  the  continuity  equation  on 
the  cross-band  component.  The  values  obtained  were  then  plotted  and  a  subjective 
isotach  analysis  was  done.  This  smoothed  the  data  by  eliminating  any  cumulus  scale 
irregularities  in  the  wind  field.  Further  smoothing  was  obtained  when  the  cross¬ 
band  speeds  were  obtained  from  the  isotach  analysis  and  applied  to  rectangular 
grid  points,  and  then  horizontal  convergences  were  computed  over  grid  spacings 
of  2.6,  5.2,  and  10.4  km.  These  divergences  were  then  integrated  upward  from 
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the  surface  in  one  third  km  steps.  A  vertical  velocity  of  zero  was  assumed  at  the 
surface.  The  vertical  velocities  obtained  were  then  combined  with  the  cross-band 
components  to  obtain  winds  that  represented  the  instantaneous  airflow  in  a  vertical 
plane  oriented  perpendicular  to  the  rainbands.  Examination  of  the  results  showed 
that  the  figures  obtained  with  a  2.6  km  grid  spacing  were  too  noisy,  and  the  ones 
obtained  with  a  10.4  km  spacing  were  too  smooth.  Thus,  Hobbs  et  al.  (1980)  used 
the  5.2  km  grid  spacing  for  analysis.  The  original  assumption  of  nondivergence  in 
the  direction  parallel  to  the  rainbands  was  checked  by  radar.  Hobbs  et  al.  (1980) 
compared  Doppler  measurements  in  the  north-northeast  and  the  south-southwest, 
that  is,  the  directions  parallel  to  the  frontal  system.  The  magnitude  of  divergences 
found  were  almost  always  less  than  1.0  X  10-4  s_1.  The  cross-band  values  ranged 
from  3  x  10-4  s_1  to  19  X  10~4  s-1.  These  results  indicate  that  the  rainbands  are 
approximately  two-dimensional,  and  that  the  general  flow  pattern  in  the  rainbands 
can  be  found  by  using  the  cross-band  velocity  component. 

The  divergence  field  of  one  of  the  wide  bands  was  examined  after  it  passed  the 
surface  front.  At  the  ground  a  very  shallow  layer  of  convergence  was  seen.  Above 
this,  up  to  3  km,  the  lower  layer  was  characterized  by  weak  divergence.  From  3-5  km 
in  height  there  was  an  area  of  convergence  with  magnitude  of  5  x  10~4  s-1.  This 
created  ascent  aloft,  with  vertical  velocities  of  approximately  1  ms-1.  The  airflow 
through  the  rainband  was  generally  from  the  front  to  the  rear.  The  updraft  area  in 
the  narrow  band  was  about  5  km  wide,  and  up  to  the  height  of  1.5  km  with  speeds 
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of  1.2  ms-1.  The  updraft  originated  in  an  area  of  strong  low-level  convergence, 
with  values  of  up  to  1.9  x  10“3  s_1.  The  air  that  flowed  parallel  to  the  front  in  the 
lowest  500  m  moved,  in  a  relative  sense,  toward  the  surface  pressure  trough,  where 
it  ascended  in  the  narrow  updraft.  When  this  updraft  reached  the  3  km  height, 
strong  upper-level  westerly  winds  moved  the  precipitation  eastward.  A  downdraft 
was  coupled  with  this  updraft.  This  downdraft  was  also  about  5  km  in  width  and 
had  speeds  of  approximately  0.5  ms"1.  It  was  located  between  the  heights  of  1.5 
and  3  km,  and  at  its  base  the  divergence  reached  a  value  of  1.0  x  10~3  s-1.  The 
air  within  the  cold  airmass  had  little  motion  relative  to  the  front,  but  the  air  above 
it  flowed  backwards  and  ascended  relative  to  the  front.  This  ascent  was  generally 
caused  by  convergence  along  the  upper  level  front.  The  greatest  values  of  vertical 
velocities  (60  cm  s"1)  were  found  over  the  wide  rainbands  to  the  rear  of  the  surface 
front. 

1.6  Storm  Movement 

When  examing  echo  motion,  the  echoes  associated  with  tropical  squall  lines 
were  omitted  by  Houze  and  Cheng  (1977),  since  the  movement  of  these  echoes  is 
faster  than,  and  unrepresentative  of,  the  general  population.  It  was  determined 
that  the  remaining  echoes  moved  at  3-5  ms-1.  The  average  direction  of  motion 
was  compared  to  the  streamlines  of  the  850  mb  and  700  mb  levels.  From  this  it 
was  determined  that  the  echoes  were  mostly  advected  by  the  850  mb  wind  flow. 
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Some  of  the  echoes  whose  area  was  less  than  500  km2  moved  in  an  erratic  fashion, 
but  when  the  echoes  of  this  size  moved  in  a  well-defined  fashion,  they  followed 
the  850  mb  streamline  winds.  Approximately  83%  of  the  isolated  echoes  examined 
were  elongated  in  the  horizontal.  Their  orientations  were  compared  to  the  GATE 
windfield.  There  was  a  tendency  for  the  echo  lines  to  be  parallel  with  the  lower 
level  winds.  The  angles  of  the  echo  lines  compared  to  the  850  mb  streamlines  was 
generally  less  than  45°.  Houze  and  Cheng  (1977)  also  observed  that  large  groups  of 
echoes  tended  to  form  into  lines.  These  lines  were  strongly  oriented  more  with  the 
low  level  winds  compared  to  the  isolated  echoes,  having  angles  with  the  streamlines 
of  ±25°.  The  results  of  Houze  and  Cheng  (1977)  from  GATE  data  (which  is  from  the 
eastern  tropical  Atlantic  Ocean)  were  analogous  to  those  determined  from  western 
tropical  Atlantic  Ocean  data. 

James  and  Browning  (1979)  stated  that  despite  the  variation  in  sizes  and  life¬ 
times,  the  convective  frontal  elements  displayed  some  common  features.  They  were 
always  oriented  clockwise  to  the  surface  cold  front  and  the  cores  all  moved  in  the 
same  direction  with  a  component  parallel  to  the  direction  of  the  low-level  jet.  The 
structure  of  the  line  showed  no  relation  to  topography.  In  many  of  the  cases  the 
fronts  passed  over  the  mountains  in  Wales  without  apparent  modification.  This 
lack  of  topographical  effect  seemed  strange  when  it  is  considered  that  the  linear 
formation  is  a  low-level  pheno.  .  non.  Another  surprising  find  was  that  despite  the 
•  ’low  nature  of  the  line  convection,  it  produced  that  intense  rain. 
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1.7  Mid-Latitude  and  Tropical  Comparison 

As  a  whole,  studies  seem  to  support  the  conclusion  that  there  are  very  few 

features  which  distinguish  tropical  squall  lines  from  those  in  the  mid-latitudes. 

The  most  significant  difference  cited  is  that  mid-latitude  squall  lines  occasionally 

are  composed  of  individually  steady  supercell-type  storms,  tropical  squall  lines  are 

not.  Several  similar  key  features  have  been  noted,  they  are: 

(1)  The  wind  field  appears  to  be  stationary  when  compared  to  the  squall  line 
reference  frame. 

(2)  There  is  a  nearly  two-dimensional  air  circulation  over  part  of  the  squall  line. 

(3)  Two  horizontal  flows  are  observed:  one  is  a  flow  toward  the  squall  line  caused 
by  tl  lines  motion  relative  to  the  environmental  motions,  and  the  other  is  a 
rear-to-front  flow  below  4  km. 

(4)  Intense  updrafts  are  located  at  the  leading  edge. 

(5)  The  downdrafts  are  weak  and  do  not  appear  to  play  a  significant  role  in  the 
system’s  dynamics. 

(6)  Leading  edge  forcing  is  caused  by  a  cold  pool  and  low-level  convergence  which 
are  maintained  by  the  rear-to-front  flow  during  the  mature  phase  of  the  squall 
line  system. 

(7)  The  leading  edge,  referiv  >  as  the  “convective  region”,  is  the  area  of  the 


maximum  reflectivity  values. 
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(8)  The  area  under  the  trailing  anvil  referred  to  as  the  “stratiform”  region  exhibits 
lower  and  more  uniform  reflectivity  values,  and  may  have  a  radar  bright  band 
(Ray,  1988). 
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Chapter  2 
Methodology 

2.1  Doppler  Radar  Analyses 

Two  main  types  of  synthesis  were  performed.  They  were:  dual-Doppler  analy¬ 
sis  and  Velocity  Azimuth  Display  (VAD).  Each  of  these  is  described  in  this  chapter. 

2.1.1  Dual-Doppler  Analysis 

The  TAMEX  data  used  in  the  dual-Doppler  analyses  was  collected  by  the 
NCAR  CP-4  and  the  NOAA  TOGA  radars.  The  wavelengths  (A)  of  these  radars  are 
5.49  cm  and  5.38  cm,  respectively.  The  beamwidth  for  the  CP-4  radar  was  1.02° and 
the  beamwidth  for  the  TOGA  radar  was  1.60°.  The  Pulse  Repetition  Frequency 
(PRF)  was  1250  Hz  for  both  radars.  The  maximum  unambiguous  Doppler  velocity, 
which  defines  the  Nyquist  co-interval,  is  given  by 

.  (i) 

Thus,  the  NCAR  CP-4  has  a  maximum  unambiguous  Doppler  velocity,  or 
Nyquist  co-interval,  of  ±  17.16  ms-1,  and  the  NOAA  TOGA  radar  has  a  Nyquist 
co-interval  of  ±  16.81  ms"1.  All  velocities  exceeding  this  interval  are  aliased 
(folded),  and  are  indistinguishable  from  the  ones  in  the  fundamental  Nyquist  co¬ 
interval,  if  conventional  Doppler  processing  is  used  (Ray,  1988).  For  example,  a 
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wind  of  20  ms-1  toward  the  NCAR  radar  would  appear  as  14.32  ms-1,  because  it 
exceeds  the  Nyquist  co-interval.  The  raw  data  contained  folded  velocities  and  noise 
that  would  create  errors  in  synthesis.  Thus,  these  data  were  manually  edited  using 
a  version  of  the  RDSS  editor  ported  to  a  Sun  UNIX  system. 

Once  all  necessary  deleting  and  unfolding  was  accomplished,  the  data  was  ana¬ 
lyzed  using  a  dual-Doppler  synthesis  method  similar  to  that  described  in  Ray  et  al. 
(1980).  First,  the  analysis  domain  was  chosen.  The  NCAR  CP-4  radar  was  chosen 
to  be  the  origin  of  the  grid  system.  For  this  system,  the  NOAA  TOGA  radar  was 
then  located  at  the  Cartesian  coordinates  of  (-14.7,  -41.3,  0.2)  km,  using  a  Lambert 
conformal  map  projection  with  23.50°N  and  25.50°N  as  reference  latitudes.  The 
grid  spacing  used  was  1.0  km  in  all  directions  and  the  spherical  influence  radius 
used  in  the  objective  analysis  was  also  1.0  km.  The  grid  domain  was  a  square, 
90  km  on  a  side  and  20  km  in  the  vertical.  Three  versions  of  this  grid  size  were 
used  for  the  dual-Doppler  synthesis.  The  location  of  these  grids  was  dependent 
upon  the  location  of  the  data.  In  each  case,  the  NCAR  CP-4  was  located  at  the 
origin.  The  grids  were:  TAMEX  West,  a  grid  designed  to  encompass  the  western 
dual-Doppler  analysis  area  over  the  Straits  of  Taiwan.  The  lower  left-hand  comer 
point  of  this  grid  is  defined  by  the  coordinates  (-70.0,  -00.0,  0.0).  TAMEX  East,  a 
grid  designed  to  cover  the  eastern  dual-Doppler  analysis  area  over  the  central  and 
eastern  part  of  island.  The  lower  left-hand  corner  point  of  this  grid  is  defined  by  the 
coordinate  point  (-20.0,  -70.0,  0.0).  TAMEX  central,  a  grid  designed  to  encompass 
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a  central  area  that  covers  parts  of  both  TAMEX  West  and  TAMEX  East.  The 
lower  left-hand  corner  point  of  this  grid  is  defined  by  the  coordinates  (-45.0,  -60.0, 
0.0). 

A  radar  collects  data  in  a  spherical  coordinate  system.  Each  sample  volume 
is  located  at  a  given  range,  azimuth,  and  elevation  from  a  radar.  The  radar’s 
beam  samples  the  atmosphere  at  regularly-spaced  intervals  along  the  beam.  These 
sample  volumes  are  called  gates,  and  are  defined  by  the  beam  width  and  a  spacing 
(typically  150-250  meters)  along  the  radial.  For  ease  of  analysis,  the  data  which  are 
collected  in  spherical  coordinates  are  converted  to  Cartesian  coordinates.  During 
the  time  that  it  takes  to  scan  through  a  storm  the  storm  evolves.  For  the  purpose 
of  analysis,  steady  state  is  assumed.  However,  the  data  are  adjusted  to  account  for 
storm  translation. 

The  reflectivity  values  are  thresholded  at  0  dBZ,  where  dBz  is  a  unitless  quan¬ 
tity  defined  by 

dBZ  =  10  x  log10  ( Jj  )  ■  (2) 

The  term  Zt  represents  the  effective  radar  reflectivity  factor  (mm6  m~3),  and  Z0  is 
1  mm6  m~3.  Power,  not  the  logarithm  of  power,  is  averaged  in  the  synthesis.  This 
conversion  is  performed  in  order  to  avoid  the  low  bias  of  the  values  that  would  result 
if  averaging  was  in  the  logarithm  of  power.  However,  all  displayed  reflectivities  are 


in  dBZ. 
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The  data  are  distance-weighted  averaged  using  the  weighting  function  given 
by  the  Cressman  (1959).  The  Cressman  weighting  function  is  mathematically  ex¬ 
pressed  as 


V/i  = 


(3) 


Here  Wi  represents  the  weight  assigned  to  an  observation  at  a  grid  point.The 
variable  d{  is  the  observation  to  grid  point  distance  and  R{  is  the  influence  radius, 
that  is,  the  maximum  distance  allowed  between  a  grid  point  and  an  observation 
point  to  which  a  weight  is  assigned.  For  this  analysis,  the  value  of  Ri  is  1.0  km. 

By  using  data  from  two  radars  and  the  continuity  equation  it  is  possible  to 
uniquely  determine  the  three  orthogonal  components  a,  v,  and  w  of  motion  (See 
Appendix  A).  These  correspond  to  the  x,  y,  and  z  directions,  respectively.  The 
equations  for  these  components  of  motion  are 


RiVliy  -  V2)  -  H2k2(y  -  yi)  -  W[(z  -  n)(y  -  y2)  -(z-  z2)(y  -  yQ]  ^ 
(x  -  xi)(y  -  y2)  ~  (x  -  x2)(y  -  yi)  ’  ’ 

RiV'iix  -  xi)  -  RjVi(x  -  xi)  -  W[(z  -  z2)(x  -  gi)  -(z  ~  zi)(x  -  x2)\  ,  . 

(x  -  x2)(y -yi)  ~  (x  -  xx)(y  -  y2)  ’  ’ 


„  du  dv  dw 

0  =  7T  +  +  n - KW 

ox  oy  oz 


(6) 


In  these  equations,  the  subscripts  denote  the  radar  number,  the  point  (jcj,  ylt 
zj)  denotes  the  location  of  the  radar  one  (NCAR),  and  (x2,  y2,  z2)  denotes  radar 
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two  (TOGA).  The  coordinate  triple  (a:,  y,  z)  is  the  grid  point  location,  k  is  the 
logarithmic  rate  of  change  of  density  with  height,  V\  and  V2  are  measured  radial 
velocities,  Hi  and  i?2  are  measures  of  range  from  each  radar  respectively,  and  W  is 
described  by  the  equation 

W  =  w  +  Vt  (7) 

Here,  Vt  is  the  particle  terminal  velocity.  We  assign  a  value  to  Vt  by  using  the 
adjusted  Joss  and  Waldvogel  (1970)  equations  (See  Appendix  A).  Then,  using  a 
value  of  zero  for  w  for  the  boundary  condition,  values  of  W  are  computed.  After 
calculating  u  and  v,  a  new  value  of  w  is  found  using  the  inelastic  continuity  equation 
(Eq.  6).  This  equation  is  integrated  down  from  the  storm  top,  where  the  initial 
boundary  condition  was  applied.  This  gives  a  new  value  of  w,  which,  through 
Eq.  7,  gives  a  new  value  of  W.  Through  Eqs.  4  and  5,  the  computation  is  iterated. 
The  iteration  continues  until  the  values  of  u,  v,  and  w  converge,  typically  just  a  few 
times.  The  process  is  repeated  at  each  level,  using  the  value  of  w  from  the  previous 
level  as  the  first  guess.  The  calculated  components  of  motion  are  then  variationally 
adjusted  to  an  imposed  constraint  that  sets  the  value  of  the  vertically-integrated 
horizontal  divergence  (See  Appendix  B). 

2.2  Errors 

Errors  are  introduced  in  the  collection  and  analysis  of  the  Doppler  velocities 
from  several  sources.  One  source  is  the  rotation  of  the  antenna.  The  rotating 
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antenna  means  that  gates  are  displaced.  This  creates  overlapping  areas  that  are 
averaged  and  one  value  is  assigned  to  a  given  gate.  This  averaged  value  may  not  be 
representative  of  the  Doppler  velocity  in  the  space  the  gate  occupies.  Interpolation 
to  a  grid  point  is  another  source  of  errors.  The  velocity  at  the  grid  point  is  an 
average  of  nearby  measured  values.  This  interpolated  value  will  vary  from  the  true 
velocity  at  the  grid  point.  Temporal  errors  can  be  caused  by  the  advection  and 
evolution  of  the  storm.  Some  advection  errors  can  be  corrected  by  inducing  storm 
motion  into  the  analysis. 

Both  radar  locations  and  the  baseline  (the  line  joining  the  radar  location)  are 
visible  in  all  wind  field  plots.  There  are  no  wind  vectors  close  to  the  baseline  due  to 
large  errors  that  arise  from  the  sampling  geometry.  The  variances  of  the  uncertainty 
in  the  wind  components  due  to  sampling  geometry  is  given  by  au2,  o2 ,  and  aw2, 
respectively  (See  Appendix  C). 


2  [Ri(y  -  y2)]Vv'12)  +  [Rijy  -  yi)}2{ov2) 

[(x  -  Xi)(y  -  y2)  -  (*  -  x2){y  -  i/i )]2 


(Kg -  *2 )(y  -  vi)  ~(z~  )(y  -  Vi )FK^2) 

((x  -  Xi)(y  -yt)-{x-  x2){y  -  2/1  )j2 


2  tggfe  -  :g2)]2(q'V-12)  +  [^l(g  -  xi  )]2(^1/12) 
((x  -  X2)(y  -  2/1 )  -  (x  -  xx )(y  -  1/2 )] 2 
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To  examine  more  closely  the  spatial  distribution  of  uncertainty,  location  of 
radar  1  is  placed  at  the  origin,  then  all  geometric  variables  (xj,  y j,  and  z\ )  axe  zero 
for  a  given  level  of  analysis.  Further,  if  radar  2  is  located  on  the  x  axis,  then  y2  and 
z2  are  also  zero.  Note  that  in  this  configuration,  the  x  axis  is  the  baseline.  This 
reduces  Eqs.  8  and  9  to 


l(*)(y)  -  (*  -  x2){y)}2 

2  _  LMf.  ~  ^)l2(gb2)  +  [•Rl(g)]a(<yVi2) 
[(x  -  x2)(y)  -  (*)(y)ja 


(11) 

(12) 


It  is  obvious  that  for  a  point  on  the  baseline  (x  axis)  the  value  of  y  is  zero  and 
the  variances  au 2  and  av2  in  Eqs.  11  and  12  are  unbounded.  Thus,  the  uncertainty 
variances  due  to  geometry  are  unbounded.  Generally,  if  the  point  measured  is  close 
to  the  baseline  (small  y  values),  the  error  variances  are  very  large.  This  large 
variance  at  points  close  to  the  baseline  is  what  makes  the  near  baseline  plotted 


vectors  suspect. 
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Equation  10  is  dependent  upon  the  values  of  the  error  variances  of  u  and  v. 
The  error  variance  of  w  is  either  unbounded  or  suspect  at  or  near  the  baseline, 
respectively. 

2.3  Velocity  Azimuth  Display 

One  method  of  using  the  measurements  of  horizontal  winds  from  a  radar  scan 
is  the  Velocity  Azimuth  Display  (VAD)  (Lhermitte  and  Atlas,  1961).  When  the 
wind  and  precipitation  fields  are  horizontally  homogeneous,  a  Doppler  radar  may 
be  used  to  measure  both  wind  particle  and  fall  velocity.  Given  a  tilted  radar  beam, 
the  radial  velocity  is  caused  by  both  the  particle  fall  velocity  and  the  motion  of  the 
particle  with  the  horizontal  wind.  Browning  and  Wexler  (1968)  used  VAD  tech¬ 
nique  to  analyze  the  harmonics  of  the  VAD.  The  mean  of  the  Doppler  velocity,  or 
the  magnitude  of  the  “zeroth”  harmonic,  permits  an  estimate  of  the  divergence. 
The  first  moment  will  be  displayed  as  a  sine  curve.  The  amplitude  of  this  curve  will 
give  a  measure  of  the  maximum  velocity  at  the  range  at  which  the  winds  are  exam¬ 
ined.  The  phase  of  the  first  harmonic  gives  the  direction  of  the  wind  at  the  sampled 
altitude.  Lhermitte  and  Atlas  (1961)  concurred  with  this  by  stating  that  the  di¬ 
rection  of  the  horizontal  wind  was  given  by  the  azimuth  of  the  maximum  velocity. 
The  displacement  of  the  entire  curve  from  the  zero  of  velocity  is  a  measurement 
of  the  precipitation  fall  speed  or  divergence.  The  second  spectral  moment  yields 
resultant  deformation  and  orientation  axis  of  dilatation.  Browning  and  Wexler’s 
VAD  analysis  was  done  for  the  cases  where  neither  the  wind  velocity  field  nor  the 
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precipitation  fail  speed  field  were  horizontally  uniform,  as  was  done  by  Lhermitte 
and  Atlas  (1961).  It  was  proposed  that,  provided  certain  precautions  were  main¬ 
tained,  the  effects  of  the  inhomogeneities  in  fall  speeds  and  vertical  wind  shear  will 
be  minimized.  The  method  of  VAD  analysis  used  follows  Browning  and  Wexler 
(1968)  and  is  found  in  Appendix  D. 

Lhermitte  and  Atlas  (1961),  Browning  and  Wexler  (1968),  and  Rabin  and  Zrnic 
(1980)  all  discuss  the  problems  involved  with  using  the  radial  velocity  values  when 
the  elevation  angles  are  greater  than  10°.  At  elevation  angles  above  this,  the  radial 
velocity  component,  given  by  the  fall  velocity,  produces  significant  error  in  the 
interpretation  that  the  radial  velocities  represent  the  horizontal  wind  flow.  At 
low  elevation  angles  this  fall  velocity  component  is  not  a  significant  factor  in  the 
measurements  of  radial  velocity;  that  is,  the  radar  return  represents  the  horizontal 
wind  component. 

The  best  way  to  avoid  these  errors  caused  by  inhomogeneities  is  to  employ 
an  optimum  scanning  procedure.  This  procedure  was  outlined  by  Browning  and 
Wexler  (1968).  They  proposed  that  for  a  rain-only  (no  frozen  precipitation  present) 
situation,  the  elevation  angle  should  not  exceed  9°,  and  that  the  horizontal  distance 
over  the  earth’s  surface  should  not  exceed  20  km.  When  examining  cases  that 
involve  nonuniform  rail  appreciable  amounts  of  convection,  these  restrictions  do 
not  apply.  The  elevation  angle  maximum  is  set  as  a  result  of  the  errors  introduced 
by  fall  speed  inhomogeneities,  and  the  range  maximum  is  set  as  a  result  of  the 
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errors  introduced  by  vertical  wind  shear.  Browning  and  Wexler  (1968)  do  state, 
however,  that  if  there  is  reason  to  believe  that  the  fields  examined  are  relatively 
homogeneous  that  the  limits  can  be  relaxed.  The  choice  of  operating  criteria,  given 
the  limitations,  is  a  compromise  among  three  considerations: 

1.  Higher  altitudes  can  be  observed  only  at  longer  slant  ranges. 

2.  The  accuracy  of  estimates  of  divergence  is  enhanced  by  making  the  zeroth 
harmonic  coefficient  as  large  as  possible,  and  this  is  done  by  using  longer  ranges. 

3.  A  large  percentage  of  the  scan’s  range  ring  examined  must  be  filled  with  useable 
velocity  values.  This  tends  to  occur  at  shorter  ranges. 

Browning  and  Wexler  (1968)  further  suggested  a  method  of  determining  the  accu¬ 
racy  of  the  technique  by  comparing  the  results  from  different  elevation  Jingles  where 
the  different  ranges  are  at  the  same  altitude. 
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Chapter  3 
Analysis 


3.1  Dual-Doppler  Analysis 

Dual-Doppler  analysis  at  082029  LST  (Figs.  9a-g)  depicts  the  leading  areas 
of  the  echo  as  it  enters  the  dual-Doppler  radar  analysis  network.  At  the  surface 
(Fig.  9a),  a  small  echo  is  observed  in  the  southwestern  corner  of  the  plot.  The  main 
body  of  the  echo  is  to  the  west  of  the  area  depicted.  The  echo  core  intensity  is 
greater  than  30  dBZ,  and  the  wind  flow  is  from  the  southwest.  All  plots  show  the 
outline  of  Taiwan,  the  NCAR  and  TOGA  radar  locations.  Only  intensities  greater 
than  20  dBZ  are  indicated.  In  some  figures  there  arc  echoes  that  do  not  have  velocity 
vectors  associated  with  them.  This  is  caused  by  a  lack  of  velocity  data  from  one  or 
both  of  the  radars.  Velocity  data  from  both  radars  is  required  for  the  dual-Doppler 
velocity  synthesis.  The  echo  at  0820  LST  (Fig.  9b)  is  not  associated  with  the  front, 
since  at  this  time  the  front  is  still  to  the  northwest  of  this  analysis  domain.  This 
echo  is  caused  by  the  confluence  between  the  mesoscale  ridging  over  the  island 
(Trier  et  al.  1989)  (Fig.  4)  and  the  synoptic  scale  flow  from  the  approaching  frontal 
trough  (Fig.  5).  Flow  remains  from  the  southwest  up  through  the  3.0  km  height 
(Figs.  9a-d),  then  it  begins  to  veer  to  a  more  zonal  flow  (Figs.  9e-g).  This  veering 
of  the  winds  with  height  is  an  indication  of  warm  air  advection.  The  flow  pattern 
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Fig.  9a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  082029,  i 
the  surface.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-han 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Rada 
positions  and  the  coastline  of  Taiwan  are  also  given. 
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^  9h'  ^[  Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  082029  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  9c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  082029,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fiq.  9c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  08S0S9,  at 
4-0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  9f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  082029,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  9g.  Dual-Doppler  analysts  of  reflectivities  and  horizontal  velocities  for  082029,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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also  corresponds  to  the  synoptic  scale  flows  seen  in  the  upper  levels  of  Fig.  5.  The 
wind  speed  is  approximately  13  ms-1.  The  echo  core  increases  in  size  from  the 
surface  to  the  2.0  km  height  (Fig.  9a-c),  then  decreases  in  area  until  it  vanishes 
between  the  5.0  km  and  6.0  km  height  (Fig.  9e-f).  The  top  of  the  echo  is  just  over 
the  6.0  km  height. 

The  dual-Doppler  analysis  at  083559  LST  (Figs.  lOa-h)  also  shows  southwest¬ 
erly  flow  at  low  heights  (Figs.  10a-d),  that  turns  westerly  aloft  (Fig.  lQe-g).  The 
storm  echo  in  the  southwest  increased  in  area  from  0820  LST,  and  the  tops  extend 
up  to  7.0  km  in  height.  At  the  surface  (Fig.  10a)  several  small  echo  cores  greater 
than  30  dBZ  axe  observed.  These  increase  in  size  until  they  merge  into  one  large 
core  at  the  2.0  km  height  (Fig.  10a~c).  It  is  at  this  height  that  the  largest  core 
area  is  observed.  The  largest  core  area  at  082029  LST  was  also  observed  at  the  2.0 
km  height  (Fig.  9c).  The  large  core  at  2.0  km  breaks  up  into  smaller  cores  at  the 
4.0  km  height  (Fig.  le)  and  tops  at  just  over  the  5.0  km  height  (Fig.  lOf).  Weak 
convergence  at  the  surface  and  1.0  km  height  is  found  below  this  core.  Divergent 
flow  is  found  above  this  core  at  the  upper  levels  (Figs.  lOf-g).  There  is  also  more 
echo  area  observed  to  the  north  of  the  southwestern  echo  area  (north  of  y  —  —33). 
During  this  time  the  mesoscale  ridge  is  undergoing  a  sharpening  in  its  configuration 
as  the  synoptic  scale  front  approaches.  This  ridge  sharpening  continues  until  the 
front  arrives  (Figs.  4,  6,  and  7).  The  echo  propagation  speed  was  calculated  to  be 
13.4  ms-1. 
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Fig.  10a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  088559  at 
the  surface.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  10b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  083559,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Ftg.  10c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  08S559,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  lOd.  Dual  -  Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  083559,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right  hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taivian  are  also  given. 


REFLECTIVITY  (DBZ) 


03 


005*02 

Mirren 


Fig.  10 f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  08S559,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  lOh.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  088559,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Figures  lla-g  depict  the  analysis  at  085133  LST.  The  front  still  has  not  entered 
the  analysis  domain.  Since  083559  LST,  the  echo  propagated  toward  the  NCAR 
radar  another  13.5  km,  giving  a  speed  of  14.3  ms-1.  The  average  horizontal  wind 
speed  in  the  echo  area  is  approximately  15  ms-1.  Several  echo  cores  greater  than 
30  dBZ  are  now  found  at  the  surface  (Fig.  11a)  in  the  southwestern  echo  area.  As 
before,  these  cores  merge  into  a  larger  core  aloft,  which  has  its  greatest  extent  at 
the  2.0  km  level  (Fig.  11c).  Above  this  height  the  core,  for  the  first  time,  stays 
relatively  unified  as  a  single  area  (Figs,  lld-f).  The  top  of  the  echo  ore  now 
reaches  6.0  km  in  height  (the  echo  tops  are  still  just  over  7.0  km  in  height).  This 
increase  in  core  height,  along  with  a  larger  core  area,  indicates  that  the  echo  core 
has  grown.  Again,  associated  with  this  echo  core  is  convergence  near  the  surface, 
(Figs,  lla-c)  and  divergence  aloft  (Figs.  1  If— h) .  The  average  winds  continue  to  veer 
(from  southwest  to  west)  with  increasing  height.  In  the  upper  levels  (Figs.  Ilf— h) 
at  the  northwestern  portion  of  the  wind  vectors  associated  with  the  southwestern 
echo  (near  y  =  —30),  the  winds  do  not  turn  westerly.  This  is  most  likely  due  to  the 
increasing  influence  of  the  synoptic  scale  frontal  trough. 

The  090524  LST  analysis  is  depicted  by  Figs.  12a-h.  This  analysis  displays  the 
same  general  flow  pattern  as  the  two  previous  analyses.  The  echo  propagated  11  km 
closer  to  the  NCAR  radar  since  085133  LST,  giving  it  a  speed  of  11.1  ms-1.  The 
average  analyzed  speeds  are  still  near  15  ms_l,  but  a  few  winds  exceed  18  ms-1. 
The  northerly  extending  line  has  intensified.  At  the  surface  (Fig.  12a),  a  large 
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Fig.  11a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  085188,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  Jib.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  085133,  at 
1.0  km.  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  tic.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  0851SS,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  lid.  Dual  Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  0851SS,  at 
S.O  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 


71 


DATE  060007  TIME  005133  4.0  KM  0.0  DEO 


o.o  M/a 


REFLECTIVITY  (DBZ) 


0^31 

mXlKM 


00^402 
M  wCTCR 


Fig.  lie.  Dual- Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  085 IS 8,  at 
j.O  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  Ilf  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  0851SS,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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jPtg.  llg.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  0851SS,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  12a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  12b.  Dual-Doppler  analysts  of  reflectivities  and  horizontal  velocities  for  090524,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  12c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524, 

2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  lid.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocitits  for  090524,  at 
3.0  km,  above  tht  ground.  The  length  scale  fur  ike  vectors  ij  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  12e.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524 > 

4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  12f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524, 

5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  ISg.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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-A'ty-  l£h.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  090524,  of 
7.0  Aim  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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echo  core  and  several  smaller  echo  cores.  For  the  first  time,  a  large  area  of  echo 
greater  than  30  dBZ  is  found.  Prior  to  this,  the  surface  maximum  intensities  have 
been  relatively  small  in  area  and  have  merged  together  at  the  2.0  km  height.  Wind 
vectors  near  the  coastline  axe  almost  from  the  south.  The  other  winds  are  principally 
from  the  southwest.  Thus,  there  is  confluent  flow  between  the  mesoscale  ridging  (as 
indicated  by  the  southerly  flow)  and  the  synoptic  scale  troughing  (as  indicated  by 
the  southwesterly  flow).  It  is  expected  that  this  area  would  be  conducive  for  further 
development.  This  southerly  flow,  however,  is  very  shallow,  for  even  at  the  1.0  km 
height  (Fig.  12b)  all  flow  is  from  the  southwest.  At  the  2.0  km  (Fig.  12c)  height, 
the  echo  cores  have  once  again  merged.  The  core  associated  with  the  northerly 
extension  appears  to  merge  with  the  core  in  the  main  echo.  From  this  height 
upwards  (Figs.  12d-f)  the  larger  core  area  once  again  breaks  into  smaller  areas. 
The  echo  depicting  the  line  shows  the  most  rapid  decrease  with  height.  The  core 
is  almost  nonexistent  at  the  4.0  km  height,  and  the  line  echo  itself  is  almost  gone 
by  the  5.0  km  height.  Thus,  at  this  time  the  line  is  relatively  intense,  but  shallow. 
Almost  all  of  the  echo  is  gone  by  the  7.0  km  height  (Figs.  12f-h).  The  highest  echo 
top  just  reaches  the  6.0  km  height.  In  Figs.  12f-h,  the  southwest  flow  and  westerly 
flow  that  was  present  at  085133  LST  (Figs.  1  If— h)  is  dominant. 

The  propagation  speed  of  the  echo  (towards  the  NCAR  radar)  at  095308  LST 
(Figs.  13a-h)  is  estimated  at  12.6  ms-1.  At  095308  LST,  several  horizontal  velocity 
values  exceed  17  ms-1.  Examination  of  the  lowest  heights  (surface  to  2.0  km, 
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Fig.  18a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095808,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  18b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095S08,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  13c.  Fual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095308,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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fty.  18d.  Dual-Doppler  analysts  of  reflectivities  and  horizontal  velocities  for  095808,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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F»$.  15c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095S08,  at 
f.O  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  ISf.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095S08 ,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  ISg.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095808,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Jfy.  1  Sh.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  095S08,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Figs.  13a-c)  fails  to  clearly  show  the  linear  formation  that  was  seen  at  090524  LST 
(Figs.  12a-f).  Instead,  several  reflectivity  maxima  are  seen  at  the  surface  (Fig.  13a). 
Up  to  this  time,  only  one  large  echo  area  has  been  present.  Two  large  cores  (greater 
than  30  dBZ)  and  several  smaller  cores  now  exist.  The  area  of  low-level  convergence 
between  the  southerly  and  southwesterly  flow  is  still  present.  By  the  1.0  km  height, 
(Fig.  13b),  the  flow  is  southwesterly  and  the  echo  area  has  increased  compared  to 
the  area  at  the  surface.  Fig.  13c  shows  the  2.0  km  height.  Consistent  with  the 
past  analysis  times,  the  cores  have  expanded  and/or  merged  to  give  the  largest  core 
area  yet  observed.  The  cores  and  reflectivities  maxima  begin  to  decrease  in  size 
and  area  between  the  3.0  km  and  5.0  km  heights  (Figs.  13d-f).  However,  at  the  5.0 
km  height  (Fig.  13f),  a  large  area  of  greater  than  30  dBZ  is  seen  over  the  location 
of  the  low  level  convergence.  The  cause  for  this,  since  it  is  only  at  this  height,  is 
that  it  corresponds  to  the  local  height  of  the  melting  level  (0°C),  creating  the  radar 
“brightband”.  At  this  height,  melting  will  cause  an  increase  in  the  reflectivity  as 
large  frozen  particles  begin  to  melt  and  develop  a  liquid  water  coat  that  makes 
them  more  reflective.  This  height  will  be  verified  by  the  examination  of  Range 
Height  Scans  (RHI)  later  in  this  chapter.  At  6.0  km  height  (Fig.  13g),  the  flow 
is  predominately  southwesterly.  At  7.0  km  height  (Fig.  13h),  the  flow  becomes 
westerly.  Flow  prior  to  this  was  westerly  at  lower  heights. 

Figures  14a-g  depict  the  analysis  at  105232  LST.  At  this  time  the  front  is  just 
to  the  north  of  the  analysis  domain.  There  are  several  areas  of  high  reflectivity  to 
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Fig.  14a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105282,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  arc  also  given. 
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Fij.  1^6.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105282,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Ft#.  I^c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105232,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  14d.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105B8&,  at 
S.O  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  He.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105282,  at 
4-0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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jPij.  /^/.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105882,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  ljg.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  105&81,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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the  west  and  southwest  of  the  NCAR  radar.  These  echoes  are  ground  clutter  from 
the  mountains  to  the  west  of  the  radar.  Echo  cores  greater  than  30  dBZ  increase 
in  size  to  a  maximum  area  at  2.0  km  (Figs.  14a-c),  then  decrease  in  size  between 
the  2.0  km  and  6.0  km  height  (Figs.  14d-g),  where  the  core  tops  are.  The  echo 
tops  appear  to  have  decreased  in  height,  as  none  are  seen  at  7.0  km  (not  shown), 
but  there  is  still  a  large  echo  area  at  6.0  km.  At  the  5.0  km  height,  an  increase 
in  reflectivity  is  due  to  the  brightband  effect.  Tne  flow  at  the  0.0  km  (Fig.  14a)  is 
from  the  south.  By  the  2.0  km  height  (Fig.  14b),  the  flow  pattern  has  once  again 
become  southwesterly.  Fig.  14c  depicts  the  2.0  km  level;  the  wind  vectors  along  the 
northern  border  of  the  domain  depict  southwesterly  flow.  If  the  vectors  indicating 
west  wind  at  the  1.0  km  height  are  indeed  caused  by  the  approaching  front,  then  the 
first  evidence  that  the  front  is  a  shallow  phenomenon  is  observed.  Southwesterly 
flow  is  predominantly  observed  throughout  the  remain  in  heights  (Figs.  14d-g). 

The  front  is  in  the  analysis  domain  at  113232  LST  (Figs.  15a-i).  At  the  surface 
(Fig.  15a)  the  front  is  depicted  by  a  heavy  dashed  line.  This  method  of  depicting  the 
frontal  boundary  will  be  used  throughout.  The  lack  of  wind  vectors  makes  frontal 
placement  difficult,  and  it  is  done  using  the  echo  pattern  depicted  and  the  PPI  scans 
from  the  radars  (not  shown).  The  front  is  approximately  5.0  km  to  the  north  of  the 
NCAR  radar.  Winds  from  the  south  (blowing  toward  the  front)  are  approximately 
13  ms-1,  and  some  winds  behind  the  front  (estimated  from  the  NCAR  Doppler  PPI 
scan)  are  as  high  as  18  ms-1.  A  large  echo  core  (greater  than  30  dBZ)  is  observed 
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Fig.  15a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11S2S2,  at 
the  ground.  The  length  scale  for  the  vectors  i3  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is  indi¬ 
cated  by  the  heavy  dashed  line. 
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Fig.  15b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11S£S£,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 
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/’iy.  15c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11S2S2,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 


104 


080087 


113232 


3.0  KM 


'•'s’S’S 


0.0  DEO 
0.0  M/3 


Him  i  »jn  >\i  n  1 1 1 1 1  u  jt  1 1  n 

/// ?•£/ 1  t  \  f J  UK  4  >  i  A{  fe/  *  1 1 

v//A  *  4  /  /  y///}/dv\ 

///f  Vt  W/rd//////?  t  /  >  r  / 

✓>////  ;k  Jj  A/U,  s  /  tv  -  ,  (  ,‘r* 

-- ////  f/  AApi /  *  >*/ r  r  **. 


/  A// 

s'SsrSS.  sSftfrSjfS/  7//// 

^  ''//s*** 

s  ~s/->.rJn  /  fJy//  »  /  sJ  I 

?*////(//(////  *  /✓//  r’f 

(  /  /  / /S 

■::sSV;//  s/f  ass/  Ass.,,  f' 
^sA-Jfy  s/sf  /  ttsfAAA^smm  rn 

.^S/^SSSS^^  //,/,£  ' 

XC252&/&A  :  f  '  i 

y.£r;..:y  • 

S’S’SlS S  S'ASySS  *.  A,  A*7  X;  I 


'  + 

CAR 


] 

’^A% ' 

':' *■:  KK't , 

«piil 


-70.0-  -61.0  -52.0  -43.0  -34.0  -25.0  -16.0  -7.0  2.0  11.0 

REFLECTIVITY  (DBZ) 


P, 300^*02 
MM  wCIDR 


0.3<X 

MMtIHLM 


jPt^.  J5d.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11SSSS,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  15e.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  118282,  at 
4-0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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F*^.  15/.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11SSSS ,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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J?*0.  15g.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11SSS2,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  15k.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11S2SS,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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15i.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  118282,  at 
8.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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along  the  frontal  boundary.  At  the  point  (-10, 12)  there  is  a  small  echo  greater  than 
40  dBZ.  This  is  the  first  time  reflectivities  this  large  were  observed.  The  30  dBZ 
cores  have  almost  disappeared  by  the  4.0  km  height  (Pigs.  15a-e),  but  again,  at 
5.0  km  (Fig.  15f)  brightband  data  distort  this  perception.  All'of  the  cores  along 
the  front  have  vanished  by  the  6.0  km  height  (Fig.  15g).  An  elongated  echo  extends 
from  the  front  northward  along  the  line  x  =  —5.0.  This  is  a  range-folded  echo.  This 
is  an  echo  that  is  at  a  distance  greater  than  the  maximum  unambiguous  range  that 
appears  as  an  echo  at  a  distance  less  than  the  maximum  unambiguous  range,  caused 
by  the  radar’s  inability  to  distinguish  between  a  return  from  the  most  recent  pulse 
and  more  distant  returns  from  a  prior  pulse.  These  echoes  will  appear  elongated 
in  the  plots  (and  on  the  radar  screen)  as  their  angular  width  is  unchanged  when 
displayed  at  the  apparent  range  (Ray,  1988).  An  echo  that  parallels  the  coastline 
is  also  seen.  This  line  of  echoes  is  associated  with  the  mesoscale  ridge  and  synoptic 
trough  confluent  line  as  analyzed  by  Trier  et  al.  (1989).  Echo  cores  of  greater  than 
40  dBZ  are  observed  along  this  line  from  the  1.0  km  height  to  the  3.0  km  height. 
Echo  cores  greater  than  30  dBZ  extend  up  above  6.0  km,  but  are  not  present  at 
the  7.0  km  height  (Fig.  15g  and  15h).  Not  only  do  the  echo  intensities  increase 
at  this  time,  but  the  tops  are  now  up  to  a  maximum  of  8.0  km  height  (Fig.  15i). 
To  the  northeast  of  this  line  of  echoes  there  is  weak  northeasterly  wind  flow.  This 
flow  becomes  stronger  and  southwesterly  by  the  2.0  km  height.  It  is  noted  that  the 
echoes  along  the  front  greater  than  40  dBZ  are  located  at  the  intersection  of  the 
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front  and  the  ridge/trough  confluent  line,  referred  to  as  the  “triple  point”.  At  the 
2.0  km  height  (Fig.  15c),  wind  vectors  are  observed  north  of  the  frontal  boundary. 
Since  these  vectors  are  uneffected  by  the  front  (as  they  are  still  depicting  southerly 
or  southwesterly  flow)  the  fact  that  the  front  is  shallow  is  supported.  At  the  upper 
levels,  southwesterly  flow  dominates,  except  in  the  area  just  to  the  west  of  the 
TOGA  radar,  where  the  ridging  is  strongest. 

Figures  16a-h  depict  the  dual-Doppler  analysis  at  115142  LST.  At  the  sur¬ 
face,  the  front  is  now  just  south  of  the  NCAR  radar.  Frontal  speed  between  the 
113232  LST  position  and  the  current  position  is  estimated  at  18  ms-1.  Winds  be¬ 
hind  the  front  (from  the  NCAR  PPI  scan)  have  a  maximum  value  of  20  ms-1. 
Further,  just  south  of  the  radar  is  the  area  of  highest  reflectivities  (greater  than 
40  dBZ).  This  core  is  located  in  the  approximate  area  of  the  triple  point.  It  is  also 
in  the  area  where  orographic  lifting  may  begin  to  influence  convection.  This  echo 
core  (greater  than  40  dBZ)  is  observed  up  to  the  2.0  km  height  (Figs.  16a-c).  The 
only  area  of  greater  than  40  dBZ  seen  at  the  2.0  km  height  (Fig.  16c)  is  at  the 
triple  point.  At  the  surface,  to  the  southeast  of  the  front,  a  small  shallow  area  of 
cyclonic  turning  is  observed.  This  area  is  associated  with  several  cores  of  greater 
than  30  dBZ.  By  the  1.0  km  height  (Fig.  16b),  the  cyclonic  circulation  appears  to 
break  down,  and  southwesterly  flow  begins  to  dominate.  At  this  time  precipitation 
is  developing  over  the  island  in  the  vicinity  of  the  TOGA  radar.  This  results  from 
the  southwesterly  flow  and  upslope  effects.  In  Fig.  16c,  south  and  southwesterly 
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Fig.  16a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11514%,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  ia  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is  indi¬ 
cated  by  the  heavy  dashed  line. 
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Fig.  16b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11514S,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 


t 


114 


DATE  080887  TIME  115142  2.0  KM  0.0  DEG 

0.0  M/3 

24.0 
18.0 

12.0 
fl.O 
0.0 
-8.0 
-12.0 
-18.0 
-24.0 
-30.0 
-38.0 
-42.0 
-48.0 

-54.0 
-80.0 

-70:0  -81.0  -52.0  -43.0  -34.0  -25.0  -16.0  -7.0  2.0  11.0 

REFLECTIVITY  (DBZ) 

0.300V02 

mxiilM  were* 

Fig.  16c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  11514&,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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fig.  lCd.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1151ft,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fij.  16c.  Dual-Doppler  analysts  of  reflectivities  and  horizontal  velocities  for  115142,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  16f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  16g.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Ftj.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  115142,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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flow  predominates.  In  the  vicinity  of  the  TOGA  radar,  the  area  of  reflectivities 
exceeding  30  dBZ  is  growing.  As  before,  the  echoes  of  this  intensity  their  maximum 
areal  extent  at  2.0  km  in  height.  The  echoes  associated  with  the  upslope  effects  do 
not,  however,  exhibit  this  characteristic.  These  echoes  continue  to  increase  in  size 
up  to  the  5.0  km  height  (Figs.  16d-f).  It  is  possible  that  the  echoes  only  increase 
up  to  the  4.0  km  height,  and  the  illusion  of  continued  size  growth  is  enhanced  by 
brightband  effects.  The  core  size  rapidly  decreases  by  the  6.0  km  height  (Fig.  16g). 
At  the  7.0  km  height,  almost  all  the  echo  has  vanished,  and  only  a  few  small  areas 
greater  than  20  dBZ  remain.  Southwesterly  wind  flow  in  the  upper  levels  is  clearly 
dominate. 

Analysis  at  123013  LST  is  shown  in  Figs.  17a-i.  At  the  surface  (Fig.  17a)  the 
front  is  clearly  observed  in  the  wind  field.  Winds  behind  the  front  (having  speeds 
of  approximately  19.0  ms-1)  are  northerly,  while  winds  ahead  of  the  front  remain 
southerly.  The  western  portion  of  the  front  is  observed  to  have  moved  southward 
(compared  to  the  position  at  115142  LST)  (Figs.  17a  and  17b)  at  approximately 
18  ms-1.  Conversely,  the  portion  of  the  front  directly  south  of  the  NCAR  radar  has 
moved  at  less  than  1.0  ms-1,  and  the  fronts  to  the  east  and  northeast  of  the  NCAR 
radar  are  virtually  stationary.  The  lack  of  movement  in  these  areas  is  caused  by 
the  shallow  front  being  blocked  by  island  terrain.  Near  the  front,  just  off  the  island 
coast,  echoes  exceed  40  dBZ.  These  cores  reach  their  maximum  area  at  1.0  km 
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Fig.  17a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  112S01S, 
at  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right- 
hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure. 
Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is 
indicated  by  the  heavy  dashed  line. 
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A'^.  17b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1SS01S,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line.  Cross-sectional  plot  given  in 
Figure  19  is  denoted  by  line  A-B. 
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F»«7.  J7c.  Dual- Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1SS01S,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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fty.  17d.  Dual-Do^ler  analysis  of  reflectivities  and  horizontal  velocities  for  123018,  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  17c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1SS01S ,  at 
4-0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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F*^.  17/.  Dual- Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  128018,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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F*0.  i7jr.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  12S01S,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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/fy.  17h.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  12S018,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  17 i.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  12S01S,  at 
8.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  arc  also  given. 


130 


height,  but  axe  observed  up  to  6.0  km  height  (Figs.  17a-g).  The  areal  extent  of  cores 
greater  than  30  dBZ  also  maximizes  at  1.0  km  height  (Fig.  17b).  These  cores  extend 
to  a  height  of  6.0  km  (Fig.  17g).  The  echoes  (greater  than  20  dBZ)  have  increased  in 
height,  with  some  tops  reaching  the  8.0  km  (Fig.  17i).  The  upsloping  echo  area  (to 
the  northwest  of  the  TOGA  radar)  is  seen  at  the  1.0  km  height.  Cores  (associated 
with  this  area  only)  greater  than  30  dBZ  are  observed  to  reach  a  maximum  areal 
extent  at  the  3.0  km  (Fig.  17d)  height.  Again,  brightband  effects  are  observed  at 
the  5.0  km  height  (Fig.  17f).  The  echoes  in  this  upslope  area  do  not  reach  heights 
greater  than  6.0  km.  It  is  also  possible  that  the  convection  in  this  area  is  not 
only  caused  by  upsloping  winds,  but  also  may  be  aided  by  a  confluence  between 
the  southerly  and  southwesterly  flow.  However,  lack  of  data  along  the  baseline 
makes  this  difficult  to  determine.  Northerly  frontal  winds  have  vanished  by  2.0  km, 
again  supporting  the  shallow  front  concept.  This  is  verified  by  a  cross-section  plot 
(Fig.  18)  corresponding  to  line  A-B  in  Fig.  17b.  In  Fig.  18,  the  frontal  flow  does 
not  reach  the  2.0  km  height  level.  Figure  18  also  shows  that  the  most  intense  echo 
core,  greater  than  40  dBZ,  is  located  near  the  leading  edge  of  the  front.  Further,  we 
can  see  from  this  cross-section  that  the  flow  is  an  anafront.  As  discussed  by  Barry 
and  Chorely  (1982),  the  frontal  area  is  characterized  by  rising  motion  in  the  warm 
sector  ahead  of,  and  over,  the  front  itself.  It  is  also  a  very  active  front,  as  indicated 
by  the  relatively  intense  echoes  associated  with  it.  The  cold  anafronts  studied  by 
Browning  and  Harrold  (1970)  had  a  narrow  band  of  convection  associated  with  a 
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broad  area  of  precipitation.  The  echo  pattern  in  Fig.  18  indicates  that  this  also 
occurred  with  this  front. 

By  124924  LST  (Figs.  19a-h),  the  western  portion  of  the  front  (relative  to  the 
analysis  grid)  has  moved  south  of  the  NCAR  radar  to  a  point  y  —  —36  (Fig.  19a). 
This  gives  fronted  speed  of  approximately  16  ms-1  for  the  western  part  of  the  front. 
The  front  lies  in  an  east  to  west  direction.  Trier  et  al.  (1989)  in  their  12  LST  analysis 
(Fig.  6)  analyzed  a  wave  on  the  front;  however,  it  is  not  seen  in  the  present  analysis. 
This  could  be  due  to  a  lack  of  data.  The  eastern  portion  of  the  front  is  oriented 
from  the  southwest  to  the  northeast.  This  part  of  the  front  is  approximately  co¬ 
located  with  the  western  edge  of  the  island  elevations  greater  than  1.0  km  (Fig.  1). 
This  gives  credibility  to  the  assumption  that  the  island  terrain  served  to  block  the 
front  as  it  moved  southward.  Further,  in  the  mesoscale  analysis  by  Trier  et  al. 
(1989),  the  frontal  areas  depicted  in  Figs.  6  and  7  also  conform  to  the  location 
of  the  island  elevations  greater  than  1.0  km.  The  upper  level  flow  patterns  at 
124924  LST  (Figs.  19d-h)  are  analogous  to  the  ones  at  123013  LST  (Figs.  17d-i). 
The  maximum  echo  tops  are  at  7.0  km  (Fig.  19h).  Only  a  few  cores  greater  than 
40  dBZ  are  seen  at  this  time,  but  none  appear  to  be  associated  with  the  moving 
portion  of  the  front.  Cores  greater  than  30  dBZ  have  tops  at  5.0  km  (Fig.  19f).  In 
the  lower  levels  (Figs.  19a  and  19b),  the  flow  is  predominately  northeasterly  behind 
the  front.  Figure  20  is  a  cross-section  through  the  front.  The  location  of  this  plot 
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Fig.  19a.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  18J9S4,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is  indi¬ 
cated  by  the  heavy  dashed  line. 
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Fig.  19b.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  ItfWt  ai 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  ia  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  w  indicated  by  the  heavy  dashed  line.  Line  A-B  depicts  the  position  of 
the  cross-sectional  plot  depicted  in  Fig.  SI. 
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Fig.  19c.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1S49&4,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 


Fig.  19d.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  lif924>  at 
3.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  19e.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1H9H,  at 
4.0  km  above  the  ground.  The  length  scale  for  the  vectors  t4  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  19f.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124924,  at 
5.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  19g.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  ttf9tf,  at 
6.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  19h.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  124918,  at 
7.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  comer.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 


Fig.  SO.  Vertical  cross-section  of  reflectivity  and  wind.  View  is  toward  the  WNW  along 
the  cross-section  (A-B)  indicated  in  Fig.  19b . 
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is  defined  by  line  A-B  in  Fig.  19b.  The  leading  edge  of  reflectivities  greater  than 
30  dBZ  is  approximately  at  the  same  location  as  the  leading  edge  of  the  front. 

At  132902  LST  (Figs.  21a-c),  the  front  is  just  north  of  the  TOGA  radar.  The 
portion  of  the  front  over  the  coast  is  moving  at  approximately  7  ms-1.  The  portion 
of  the  front  over  the  island  has  not  moved.  Near  the  frontal  boundary,  large  cores 
greater  than  30  dBZ  are  seen.  Smaller  cores  with  reflectivities  greater  than  40  dBZ 
are  seen  at  the  surface  (Fig.  21a).  These  40  dBZ  cores  have  tops  less  than  2.0  km  in 
height.  Some  of  the  30  dBZ  cores  extend  to  6.0  km  in  height,  and  some  extend  to  7.0 
km.  The  1.0  km  height  (Fig.  21b)  shows  the  orientation  of  the  cross-section  plot 
(lines  A-B)  depicted  by  Fig.  22.  The  most  interesting  feature  in  Fig.  22  is  found  at 
the  frontal  boundary  (near  the  lower  left-hand  corner  of  the  plot).  At  the  edge  of 
the  front,  winds  are  observed  to  point  toward  the  surface.  This  flow  is  co-located 
with  an  echo  core  of  greater  than  30  dBZ. 'Thus,  this  flow  most  likely  represents 
the  outflow  from  the  convective  core.  Figure  23  shows  the  1.0  km  height  frontal 
position  at  140905.  Lack  of  data  make  the  other  levels  unusable.  It  is  observed 
in  this  figure  that  the  front  still  has  not  changed  position.  Since  this  position 
is  approximately  the  same  as  the  island  elevations  greater  than  1.0  km,  terrain 
blocking  is  further  supported.  When  the  front  was  over  land,  prior  to  reaching  the 
island  elevations  greater  than  1.0  km  (where  blocking  occurred),  frontal  speed  is 
estimated  at  6  ms-1.  This  speed  corresponds  to  the  speeds  of  3-5  ms-1  found  in 
the  GATE  study  by  Houze  and  Cheng  (1977).  The  over  water  speeds  of  the 
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Fig.  21a.  Dual- Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  132902,  at 
the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower  right-hand 
corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the  figure.  Radar 
positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal  position  is  indi¬ 
cated  by  the  heavy  dashed  line. 
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Fig.  tlb.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  lSt90t,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line.  Line  A-B  denotes  the  position 
of  the  cross-sectional  plot  in  Fig.  SB. 
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Fij.  JBJc.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  182902,  at 
2.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given. 
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Fig.  tS.  Dual-Doppler  analysis  of  reflectivities  and  horizontal  velocities  for  1J0905,  at 
1.0  km  above  the  ground.  The  length  scale  for  the  vectors  is  given  in  the  lower 
right-hand  corner.  The  reflectivity  gray  scale  is  given  on  the  right  side  of  the 
figure.  Radar  positions  and  the  coastline  of  Taiwan  are  also  given.  Frontal 
position  is  indicated  by  the  heavy  dashed  line. 
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TAMEX  front  were  as  high  as  19  ms-1.  Studies  in  the  mid-latitudes  by  Browning 
and  Harrold  (1970)  and  Hobbs  ti  al.  (1980)  showed  most  of  the  fronts  to  move  in 
excess  of  10  ms-1. 

3.2  RHI  Analysis 

Figures  24-26  are  Range  Height  Indicator  (RHI)  scans  from  the  NCAR  CP—4 
radar  taken  before  (Figs.  24  and  25,  at  1100  LST)  and  after  frontal  passage  (Fig.  26, 
at  1258  LST).  In  all  of  these  figures,  the  reflectivity  field  is  at  the  top  and  the  velocity 
field  is  at  the  bottom.  The  shade  scales  that  indicate  the  values  in  the  figures  are 
located  below  each  of  the  reflectivity  and  velocity  plots.  The  range  rings  me  at 
every  10  km. 

The  depth  of  the  front  is  1.57  km.  Figure  24  is  taken  along  the  10° azimuth. 
Toward  the  left  side  of  the  reflectivity  plot,  reflectivities  of  36  dBZ  are  observed  at 
4.75  km  height.  This  is  called  the  brightband,  and  it  is  an  indicator  of  the  melting 
(0°)  level.  The  velocity  pattern  clearly  depicts  the  front.  The  leading  edge  of  the 
front  is  16.5  km  from  the  radar,  as  indicated  by  the  negative  velocity  values.  Note 
the  radar’s  location  is  at  the  left  edge  of  the  plots,  and  negative  velocities  represent 
flow  towards  the  radar,  while  positive  values  denote  flow  away  from  the  radar.  The 
front,  and  the  corresponding  cold  air  associated  with  it,  are  seen  as  a  shallow  wedge 
of  negative  velocities  underneath  the  positive  velocity  pattern  in  the  center  and 
right  side  of  the  figure.  The  maximum  radial  velocity  in  this  cold  wedge  is  15  ms-1. 
Above  the  front,  the  maximum  radial  velocity  is  16  ms-1,  at  4.3  km 


Fig.  24.  RHI  scan  of  reflectivity  and  radial  velocity  from  the  NCAR  CP-4,  (a)  is  re¬ 
flectivity  and  (b)  is  radial  velocity,  both  at  10  degrees  azimuth  at  1100S9  LST. 
Scales  are  given  at  the  bottom  of  each  figure. 
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Fig.  25.  RHI  scan  of  reflectivity  and  radial  velocity  from  the  NCAR  CP-4-  (a)  u  reflec¬ 
tivity  and  (b)  radial  velocity  for  SO  degrees  azimuth  and  110108  LST.  Range 
marks  on  the  velocity  display  btt  in  increments  of  10  km.  Scales  are  given  at 
the  bottom  of  each  figure. 


151 


height.  At  the  leading  edge  of  the  front,  the  reflectivity  pattern  indicates  that  a 
precipitation  appears  to  be  forming,  as  reflectivities  greater  than  30  dBZ  are  seen. 
A  few  kilometers  behind  the  leading  edge  of  the  front,  a  well-developed  core  of 
34  dBZ  is  seen  at  the  1.8  km  height.  To  the  right  (toward  the  rear  of  the  front) 
the  echo  weakens.  This  agrees  with  the  findings  of  Houz,  ''1977)  who  stated  that 
new  cells  form  at  the  leading  edge  of  a  line  and  mature  cells  move  toward  the  rear 
of  the  echo  pattern,  where  they  dissipate.  Between  20  and  30  km  range,  the  echoes 
(greater  than  20  dBZ)  appear  to  curve  toward  the  radar  with  decreasing  height. 
This  is  caused  by  precipitation  falling  through  a  flow  regime  that  is  initially  away 
from  the  radar  into  a  flow  regime  towards  the  radar. 


Figure  25  is  taken  along  the  30°  azimuth  at  1100  LST.  The  reflectivity  bright- 
band  is  clearly  evident.  There  are  several  cores  (32  dBZ  or  greater)  behind  the 
front.  At  lower  heights,  these  cores  merge  into  one  echo.  It  is  the  upward  extending 
sections  of  greater  than  32  dBZ,  separated  by  reflectivities  of  less  than  32  dBZ, 
that  indicate  that  this  is  a  multicelled  system.  Again,  the  cells  appear  to  form  at 
the  leading  edge  of  the  front,  increase  in  intensity  just  behind  the  leading  edge  of 
the  front,  and  then  weaken  as  they  move  toward  the  rear.  In  Fig.  25,  the  front  is 
17.5  km  from  the  radar,  and  as  in  Fig.  24  the  maximum  velocities  seen  are  15  ms"1 
toward  the  radar,  and  16  ms-1  away  from  the  radar.  In  both  Figs.  24  and  25  the 
maximum  velocity  are  also  behind  the  most  intense  and/or  largest  reflectivity  cores. 
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The  RHI  scan  at  1258  LST  is  depicted  in  Fig.  26.  By  this  time  the  front  has 
moved  south  of  the  NCAR  CP-4  radar,  and  the  orientation  (at  210°)  of  the  figure 
is  reversed  from  the  previous  RHI  scans.  The  view  is  from  behind  the  front  toward 
the  leading  edge,  which  is  almost  40  km  from  the  radar.  In  this  figure,  small  cells 
form  at  the  leading  edge  of  the  front,  and  grow  toward  the  rear  of  the  front.  It  is 
also  observed  in  Fig.  26  that  the  precipitation  trajectories  curve  as  they  approach 
the  ground.  The  maximum  horizontal  wind  speed  toward  the  radar  16  ms-1  is  the 
same  the  maximum  behind  the  front  (flow  away  from  the  radar).  The  very  leading 
edge  of  the  front  is  nearly  vertical  in  Fig.  26.  Given  this  time  and  figure  orientation, 
this  blunt  front  coincides  with  island  elevations  greater  than  1.0  km  in  height.  In 
Figs.  24  and  25,  the  front  is  over  water  or  island  elevations  less  than  1.0  km  in 
height.  In  these  figures,  the  leading  edge  of  the  front  is  slightly  tilted.  Whenever 
the  front  was  over  water  or  island  terrain  elevations  of  less  than  1.0  km  the  leading 
edge  of  the  front  had  a  slight  slope.  When  the  front’s  leading  edge  was  in  an  area 
where  the  elevations  exceeded  1.0  km  it  appeared  nearly  vertical. 

3.3  VAD  Analysis 

The  Velocity  Azimuth  Display  (VAD)  analysis  used  data  from  the  NCAR  CP-4 
radar  for  times  from  1052  LST  to  1329  LST.  All  analyses  used  an  antenna  elevation 
of  1.8°.  Slant  ranges  of  5.0  km,  10.1  km,  19.8  km,  and  29.9  km  were  examined, 
although  not  all  of  these  ranges  were  examined  at  each  time.  The  VAD  circles  were 
designed  so  that  the  VAD  circle  did  not  cut  through  the  front.  Products 
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Fig.  26.  RHI  scans  of  reflectivity  and  radial  velocity.  Radar  is  NCAR  CP-4.  (a)  is 

aUd  ^  **  Tadial  vtlociiy>  botl1  at  BIO  degrees  azimuth  at  1258 SO 
LST.  Scales  are  given  at  the  bottom  of  each  figure.  Range  marks  on  the  velocity 
display  are  in  increments  of  10  km. 
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from  this  type  of  analysis  include  wind  direction  and  speed,  convergence,  stretching 
(dilatation)  and  shearing  deformation.  Prior  to  discussion  of  the  VAD  results,  a 
brief  overview  of  these  deformations  is  given. 

In  the  lower  atmosphere,  values  of  dilatation  and  wind  shear  are  typically  of 
the  order  10~4-10~5  s'"1 .  Larger  values  are  often  found  with  strong  winds  and/or 
smaller  scales.  Correspondingly,  larger  scales  and/or  light  winds  may  decrease  typ¬ 
ical  values  by  an  order  of  magnitude.  A  stretching  deformation  field  has  two  char¬ 
acteristic  axes,  one  toward  which  the  streamlines  converge  asymptotically  (the  axis 
of  dilatation)  and  one  from  which  the  streamlines  diverge  asymptotically  (the  axis 
of  contraction).  Positive  values  of  deformation  give  the  X  axis  as  the  axis  of  di¬ 
latation,  and  Y  as  the  axis  of  contraction.  If  the  deformation  value  is  negative  the 
axes  are  reversed.  In  a  pure  deformation  field  a  parcel  of  air  is  stretched  along 
the  axis  of  dilatation  and  compressed  along  the  axis  of  contraction.  The  shearing 
deformation  field  resembles  the  stretching  deformation  field,  except  that  the  axes  of 
shearing  deformation  are  rotated  with  respect  to  those  of  stretching  deformation. 
If  the  shearing  value  is  positive  then  shear  axes  are  45°  to  the  left  of  the  axes  of 
stretching  deformation.  If  the  shear  value  is  negative  the  axes  are  45° to  the  right. 
The  air  parcels  are  distorted  through  the  wind  field  differences  in  shear.  These 
two  types  of  deformation  are  combined  into  the  total  deformation,  which  is  located 
around  one  set  of  deformation  axes  (See  Appendix  D).  Deformation  is  a  primary 
factor  in  frontogenesis.  To  a  first  approximation,  frontogenesis  occurs  when  the 
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angle  between  the  axis  of  dilatation  and  potential  isotherms  is  less  than  45°.  If  the 
angle  is  greater  than  45°,  then  frontolysis  occurs. 

The  VAD  analysis  at  1052  LST  is  given  by  Figs.  27a-d.  The  VAD  observations 
axe  given  by  Fig.  27a.  In  all  of  the  VAD  plots,  observed  velocities  are  represented  by 
crosses,  the  mean  and  first  harmonic  are  drawn  with  a  heavy  line,  and  the  thin  line 
represents  the  mean  plus  the  first  and  second  harmonics.  The  value  of  deformation 
is  3.7  x  10_4s_1  at  an  angle  of  -56°.  For  all  times,  the  deformation  was  the  greatest 
for  the  5.0  km  distance.  This  is  partly  due  to  the  effects  of  the  scale,  as  mentioned 
previously.  Data  from  the  10.1  km  VAD  circle  (Fig.  27b)  give  a  reduced  deformation 
of  3.9  x  10~5s-1.  When  the  VAD  circle  is  expanded  to  19.8  km  (note:  all  VAD 
distances  are  slant  range  measurements)  the  deformation  increases  (compared  to 
10.1  km)  to  a  value  of  1.2  x  10~5s-1.  At  this  time,  weak  divergence  (0.3  x  10-4  « 
1.7  x  10~4s_1)  is  found  in  the  5.0  km,  10.1  km  and  19.8  km  VAD  circles,  and  weak 
convergence  (1.5  x  10-5s-1 )  is  found  in  the  29.9  km  VAD  circle.  40  minutes  later,  at 
1132  LST,  the  wind  field  as  seen  by  the  VAD  analysis  changed  drastically.  The  wind 
shifted  nearly  180°.  VAD  circles  show  strong  convergence  (2.4  x  10”4  »10~3s_1) 
and  deformation  (8.7  x  10-4  wl.8  x  10-3s-1).  All  these  indicate  that  the  front 
had  just  passed  the  radar  site. 
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Fig.  27a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP~4  at 
1052.  Elevation  angle  is  1.8 Scales  are  given  at  the  bottom  of  each  figure. 
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Fig.  27b.  VAD  at  1052.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses , 
first  harmonic  by  the  'eavy  line ,  and  the  first  and  second  harmonies  together 
are  given  by  the  thin  line. 


AZIMUTH 

Fig.  27c.  VAD  at  1052.  Range  is  10.1  km.  Raw  velocity  valves  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  27  d.  VAD  at  1052.  Range  is  19.8  km.  Raw  velocity  valves  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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In  Fig.  28a,  the  radar  image  for  1152  LST  is  given.  411  VAD  circles  still  show 
convergence,  but  the  magnitude  is  smaller  than  20  minutes  ago  (2.6  x  10-4  «  4.2 
x  10-4s-1).  Deformation  is  also  weaker  (8.2  x  10~4  «»1.2  x  10-3s-1). 

The  radar  image  and  VAD  analysis  at  1230  LST  are  depicted  in  Figs.  29a-d. 
The  VAD  circles  were  done  at  5.0  km,  10.1  km,  19.8  km  (Figs.  29b,  29c,  and  29d, 
respectively).  Similar  wind  fields  were  found  at  1249  LST,  which  is  depicted  by 
Fig.  30a  (the  radar  image),  and  Figs.  30b-d  (the  VAD  circles).  The  VAD  circles  at 
5.0,  10.1  and  19.8  km  still  show  convergence,  but  the  VAD  circle  at  29.9  km  now 
shows  very  weak  divergence  (6.0  x  10~8s-1)  which  would  grow  to  2.5  x  10~6  at 
1329  LST.  It  was  observed  that  in  these  VAD  circles,  the  values  obtained,  with  the 
given  scale,  are  analogous  to  those  found  by  Saucier  (1965). 


PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-4 
1152.  Elevation  angle  is  1.8 0 .  Scales  are  given  at  the  bottom  of  each  figure. 
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Fig.  28b.  VAD  at  1152.  Range  is  5.0  km.  Raw  velocity  valves  are  given  by  ihe  crosses, 
first  harmonic  by  the  heavy  line ,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 


Fig.  29a.  PPI  scans  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP~4  at 
1230.  Elevation  angle  is  1.8° .  Scales  are  given  at  the  bottom  of  each  figure. 
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first  harmonic  hy  the  heavy  line.,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  29c.  VAD  at  1230.  Range  is  10.1  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  29  i.  VAD  at  1230.  Range  is  19.8  bn.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  30b.  VAD  at  1249.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line-,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  SOd.  VAD  at  1249.  Range  is  19.8  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmoTiics  together 
are  given  by  the  thin  line. 
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Fig.  Sla.  Wicaru  of  reflectivity  (left)  and  velocity  (right)  from  the  NCAR  CP-1  at 
1S29.  Elevation  angle  is  l.r.  Scales  are  given  at  the  bottom  of  each  figul 
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Fig.  Sib.  VAD  at  1329.  Range  is  5.0  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy' line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line . 
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Fig.  Sic.  VAD  at  1329.  Range  is  10.1  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  Sid.  VAD  at  1329.  Range  is  19.8  km.  Rato  velocity  values  are  given  by  the  crosses , 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 
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Fig.  Sle.  VAD  at  1S£9.  Range  is  £9.9  km.  Raw  velocity  values  are  given  by  the  crosses, 
first  harmonic  by  the  heavy  line,  and  the  first  and  second  harmonics  together 
are  given  by  the  thin  line. 


Chapter  4 


Conclusions 


4.1  Conclusions 

This  study  focuses  on  a  Doppler  radar  analysis  of  a  June  8,  1986  frontal  system 
from  Intensive  Operating  Period  (IOP)  #8  of  the  Taiwan  Area  Mesoscale  Experi¬ 
ment  (TAMEX).  Data  collected  from  two  Doppler  radars  was  used  for  the  analy¬ 
sis.  One  radax  was  the  National  Center  for  Atmospheric  Research  (NCAR)  CP-4 
radar,  located  at  a  site  known  as  Waipu  Beach.  The  other  Doppler  radar,  the 
National  Oceanic  Atmospheric  Administration  (NOAA)  TOGA  radar,  was  located 
at  the  CCK  Air  Force  Base,  just  west  of  Taichung,  approximately  45  km  south 
of  the  NCAR  radar.  Analysis  was  done  from  0820  LST  to  1409  LST  as  the  front 
approached  and  passed  through  the  radar  network.  Two  types  of  synthesis  were 
performed  with  the  data  set.  Dual-Doppler  analysis  was  preformed  extending  over 
an  eight  hour  period,  while  Velocity  Azimuth  Displays  (VAD)  (single  Doppler  anal¬ 
ysis)  was  preformed  for  times  when  the  front  was  in  the  vicinity  of  the  NCAR  radar. 
Also  several  vertical  cross-sections  through  the  front  were  examined. 

The  first  set  of  radar  echoes  entered  the  analysis  domain  at  0820  LST.  These 
echoes  were  observed  in  the  southwestern  portion  of  the  analysis  domain  and  they 
moved  toward  the  northeast  at  approximately  12  ins-1,  driven  by  the  synoptic  scale 
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southwesterly  winds  and  the  sharpening  mesoscale  ridge  as  analyzed  by  Trier  et  al. 
(1989).  At  1052  LST,  the  front  was  north  of  the  analysis  domain.  The  front  was  first 
observed  in  the  1132  LST  analysis.  It  was  diagnosed  as  an  anafront.  It  propagated 
south,  moving  more  rapidly  over  water  (up  to  approximately  12  ms-1)  than  over 
land.  Further,  when  the  front  moved  over  the  island,  it  became  virtually  stationary. 
This  position  of  the  front  was  approximately  co-located  with  the  area  of  the  island 
where  the  elevations  exceeded  1.0  km.  Thus,  it  appears  that  the  higher  island 
elevations  served  to  block  the  front.  This  concept  led  to  further  examination  of  the 
front.  In  the  dual-Doppler  analyses,  the  front  was  only  distinguishable  for  levels 
below  1.0  km  level.  This  was  verified  in  both  the  cross-sections  of  the  dual-Doppler 
analysis  and  in  the  RHI  scans  taken  through  the  front.  Vertical  depth  of  the  front 
was  estimated  at  1.57  km.  During  the  time  the  front  traveled  through  the  radar 
grid,  the  echo  pattern  occasionally  showed  a  “T”  pattern  in  which  the  front  would 
be  the  cross  of  the  “T”.  The  stem  of  the  “T”  was  formed  by  confluence  between 
the  mesoscale  ridging  and  the  synoptic  scale  trough.  An  intense  echo  was  observed 
at  the  intersection  of  the  “T”,  at  what  was  referred  to  as  the  “triple  point”.  This 
was  at  the  intersection  of  the  front  and  the  ridge-trough  confluence  line. 

Examination  of  the  echoes  and  their  cores  revealed  that  several  cells  were  always 
present.  Since  these  cells  were  not  always  in  the  same  place,  it  is  suggessted  that 
the  cells  are  relatively  intense  and  have  short  lifetimes.  Over  the  low  elevations  and 
the  water,  the  cells  were  observed  to  increase  in  size  with  height.  The  maximum 
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areal  extent  was  reached  at  the  2.0  km  height.  Above  this,  the  cores  decreased  in 
size,  with  the  noted  exception  of  an  occasional  abrupt  increase  in  expanse  at  the 
5.0  km  height.  This  was  caused  by  radar  brightbanding  effects,  and  verified  by  the 
NCAR  RHI  scans,  where  the  the  freezing  level  was  approximated  at  the  4.75  km 
height.  For  cells  over  the  elevations  greater  than  1.0  km,  the  echoes  followed  the 
same  height  to  size  pattern,  except  these  echoes  reached  maximum  expanse  at  the 
4.0  km  level.  Thus,  it  is  observed  that  the  echoes  reached  their  maximum  size 
approximately  2.0  km  above  ground  level.  On  the  average,  the  echoes  were  most 
often  observed  to  reach  maximum  height  between  the  6.0  km  and  7.0  km  height. 

VAD  analysis  revealed  that  prior  to  frontal  passage,  deformation  favorable  to 
frontogenesis  was  observed.  Further,  the  deformation  after  the  front  passage  was 
of  a  frontolytic  condition.  The  winds  ahead  of  the  front  were  southerly,  less  then 
10.0  ms-1  and  veered  with  height,  indicating  warm  air  advection.  Winds  behind 
the  front  were  northeasterly  and  had  an  approximate  speed  of  17  ms-1.  Prior  to 
frontal  passage,  weak  convergence  was  observed.  After  the  front  passes,  the  VAD 
circles  revealed  general  divergence. 

It  was  observed  that  the  echoes  associated  with  the  front  formed  at  the  leading 
edge  and  moved  towards  the  rear  as  they  matured.  This  agrees  with  previous  studies 
in  both  the  tropics  and  in  the  mid-latitudes. 
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Appendix  A 

Dual— Doppler  Radar  Synthesis 

The  synthesis  uses  two  radars  in  the  solution.  The  dependent  variables  desired 
are  the  components  of  the  wind  motion,  u,  v,  and  to.  The  equation  that  describes 
the  radar  observed  velocity  is: 

u(x  -  Xi)  +  v(y  -  yi) +  W(z  -  Zi)  =  R,Vi  for  *  =  1,2  .  (Al) 

Here  the  subscript  i  denotes  the  number  of  the  radars.  The  point  (xt>  y;,  z*) 
denotes  the  location  of  the  i<A  radar  and  the  point  ( x ,  y,  z)  is  the  grid  point  location. 
For  a  given  point,  Vi  is  the  measured  radial  velocity.  R,  is  a  measure  of  the  range 
from  the  radar  described  by 

i. 

Ri  -  (x-  x,)2  +  (y  -  y,)2  +  {z-  zi )2  .  (A2) 

In  Eq.  Al,  IV  is  the  hydrometeor  fall  speed  containing  the  vertical  velocity,  to, 
and  the  particle  terminal  fall  speed,  V*.  Thus,  W  is  described  by  the  equation 

W  =  w  +  Vt  .  (A3) 

At  this  point  the  unknown  variables  are  u,  v,  to,  W  and  Vt.  To  reduce  the 
number  of  unknown  variables  to  obtain  a  unique  solution  to  these  equations,  a  value 
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for  Vt  is  found.  This  terminal  fall  speed  is  found  using  the  Joss  and  Waldvogel 
(1970)  and  the  Foote  and  du  Toit  (1969)  relationships.  Foote  and  du  Toit  (1969) 
determined  fall  speeds  using  density  of  the  atmosphere.  As  density  increases  closer 
to  the  ground  the  fall  speed  of  hydrometer  decreases.  Joss  and  Waldvogel  (1970) 
used  reflectivity  to  find  the  terminal  velocity.  The  higher  the  value  of  reflectivity, 
the  higher  the  value  of  the  terminal  velocity.  This  relation  is  modified  by  the  Foote 
and  du  Tbit  (1969)  density  relationship.  The  relationship  between  terminal  velocity 
below  the  4  km  msl  levels  and  reflectivity  is  expressed 


Vt  =  -2.6(Ze)107e(-4K2)  ,  (4) 

where  Ze  is  the  effective  reflectivity  factor,  k  is  the  negative  of  logarithmic  rate  of 
change  of  density  with  height  km-1,  z  is  the  height  above  mean  sea  level  (in  km) 
at  the  grid  point.  Between  4  and  10  km,  the  relationship  for  Vt  becomes 


Vt  =  -2.6 (Ze)107e(  ^ 


(!Q  —  z)  (z-  4) 


(A5) 


Above  the  10  km  level  a  value  of  -2  ms-1  is  used  (See  Figure  Al). 

The  boundary  condition  value  of  zero  for  w  is  specified  at  the  storm  top.  The 
value  for  W  can  be  calculated  using  Eq.  A3.  Equation  1  is  now  solved  for  the 
remaining  unknowns  u  and  v,  in  terms  of  W.  At  radar  #1,  i  =  1,  and  Eq.  1 


becomes 
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TERMINAL  VELOCITY  (M/S) 


Fig.  At.  Terminal  velocities  from  the  adjusted  Joss  and  Wuldvogel  (1970)  equations. 

Plotted  lines  are  done  for:  (A)  10  dBZ,  (B)  BO  dBZ,  (C)  SO  dBs,  (D)  40  dBZ, 
and  (E)  50  dBZ.  Negative  velocities  denote  downward  speed.  All  plots  have  a 
value  of  -B  ms~x  above  the  10  km  height  level. 
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u(x  -  xx)  +  v(y  -  yx)  +  W(z  -  zx)  =  RXVX  .  (.46) 

Correspondingly,  for  radar  #2,  Eq.  1  is 

u(x -x2)  +  v(y-y2)  +  W(z  -  z2)  =  R2V2  .  (A7) 

Now  Eqs.  A5  and  A6  are  solved  for  u  by  multiplying  Eq.  5  by  (y  —  y2)  and 
Eq.  6  by  (y  —  yi).  This  gives 

u(x~x1)(y-y2)  +  v(y-y1)(y-y2)  +  W(z-z1)(y-y2)  =  RlV1(y-y2)  ,  (A  7) 
and 

u(x -x2)(y-yi)  +  v(y-y2)(y-y1)  +  W(z -z2)(y-y1)  =  R2V2(y-y1)  .  (A 8) 

Then  Eq.  A8  is  subtracted  from  Eq.  A7,  eliminating  the  v. 


u[(x  -  xj)(y  -  y2)  -(x-  x2)(y  -  yi )]  =R\Vx(y  -  y2)  -  R2V2(y  -  y i) 

-  W[(z  -  zx)(y  -  y2)  -(z-  z2)(y  -  yi )] 

(A9) 

Thus, 

u  _  R\V\(y  -  V2 )  -  RtV2{y  -  yi)  -  w[(z  -  zx){y  -  y2)  -  (z  -  z2)(y  -  yx )] 

(x  -  xi)(y  -  y2)  -  (x  -  x3)(y  -  yx) 


(A10) 
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Correspondingly,  an  equation  for  v  can  be  derived,  and  is  given  by 


R*Ux  -  Xj)  -  R\Vi(x  -  sx)  -  W[(z  -  z2)(x  -  xx)  -  (z  -  Zi)(x  -  xa)] 
(x  -  x2)(y  -  vi)  -  (*  -  an)(y  -  y2) 


(All) 


Given  a  value  for  W,  from  Eq.  A3,  the  values  of  u  and  t;  are  found  using 
Eqs.  A10  and  All.  Next,  a  new  vertical  velocity  iv  can  be  found  from  the  anelastic 
continuity  equation 


du  dv  dw 


o 


(412) 


The  newly  determined  value  of  to  is  put  back  into  Eq.  A3,  and  a  new  value  of  W 
is  determined.  The  iterative  process  continues  until  convergence  of  the  solution  is 
obtained.  After  convergence  has  been  obtained,  the  value  of  to  in  the  layer  calculated 
is  used  as  a  boundary  condition  for  calculating  the  next  lowest  layer.  This  process 
is  repeated  from  the  storm  top  to  the  surface  in  order  to  evaluate  the  entire  storm 


volume. 
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Appendix  B 

Variationally  adjusted  values  of  u,  v,  and  tv  are  frequently  used  to  provide  a 
more  accurate  representation  of  the  true  wind  field.  This  is  done  using  the  integrated 
horizontal  divergence  as  a  strong  constraint.  This  means  that  the  difference  in  the 
upper  and  lower  boundary  condition  is  specified  (Ray  et  al .,  1980).  The  deviations 
between  the  new  values  and  the  observed  values  of  the  horizontal  winds  will  be  a 
minimum.  If  these  deviations  are  normally  distributed,  a  functional  can  be  written 
as 


E[u,v,w]=J  j  a2  —  u°j  +/?2^u  —  v°^j 

+  {x[S  ,{%^%)dz+c]ixiy 


(Bl) 


The  variables  u  and  v  are  the  adjusted  components  of  motion;  the  variables 
u°  and  v°  are  the  observed  horizontal  components  of  motion.  The  variables  a2 
and  /?2  are  weights  which  are  inversely  proportional  to  the  error  variances  in  u  and 
v  respectively.  Thus,  a2  =  .  Error  variance  is  a  measure  of  the  differences 

between  the  estimates  and  the  true  values  of  u  and  v.  Clearly,  a  large  variance  value 
will  yield  a  small  weight  value.  In  the  intergrand’s  last  term,  which  is  the  vertically 
integrated  horizontal  divergence,  A  is  the  Lagrange  multiplier.  It  is  the  value  of  A 
that  we  will  determine  in  order  to  variationally  adjust  u,  v,  and  to  compute  w. 
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The  integrand  of  the  functional,  /,  is  next  used  in  a  general  set  of  Euler- 
Lagrange  equations.  The  set  of  Euler- Lagrange  equations  yield  three  equations 
and  three  unknowns.  These  unknowns  will  be  A,  u,  and  v.  This  gives  us  a  set  of 
Euler-Lagrange  equations  given  by 


d  dl  dl_ 
dx  q( +  du 

d  dl  <H_  _ 
dy  J  +  dv 


w 


d 


(*) 


dl 

dx 


=  0 


These  equations  are  solved,  giving  respectively 


d\ 


+  2aa  (u  -  u°)  ^  0 

_g+2^(v_t;O)=0 


(B2) 


(B3) 


(B4) 


(B5) 

(B6) 


/**•  /5u  di A 

l  +  0 

Next,  B5  and  B6  are  solved  for  u,  and  v  giving 


u  =  u°  + 


1  <9A 
2a2  9z 


(B7) 


m 


and 


v  =  u  + 


1  dX 
2pdy 


{B 11) 


/ 
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These  values  of  u,  u,  can  now  be  substituted  into  Eq.  B7,  yielding  an  elliptic 
partial  differenced  equation  which  now  has  only  one  variable,  A. 


nd_  /  1  <9A\  d_  f  1  d\\ 
dx  \2a2  dx)  +  dy  \2/32  dy  J 


rZT  (du«  dv°\  J  „ 

=  -Jo  (a7  +  w)*  +  c 


(B 12) 


FVom  Eq.  B12  a  value  of  A  is  determined.  This  can  then  be  used  in  Eqs.  BIO, 
Bll,  and  B12,  to  get  the  adjusted  values  for  u,  v,.  These  adjusted  values  will  now 


satisfy  the  divergence  equation. 
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Appendix  C 


Uncertainty  in  Wind  Field  Estimates 


FVom  Appendix  A,  the  u,  v,  and  w  components  of  the  wind  motion  are  given 
by  the  equations 


_  R\V\(v  -  V2 )  -  R2Vi{y  -  V\)  -  W[(z  -  zi)(y  -  y2)  -  (*  -  z*){y  -  yi)l  (ru 
(x  -  Xi)(y  -  y2)  -  (x  -  x2)(y  -  yi)  ’  1 

RiVi(x  -  Ji)  ~  RiVijx  -  xQ  -  W[(z  -  z2)(x  -  Xj)  -(z-  z\)(x  -  x2)]  (ro, 
V  (x  -  x2)(y  -  yx)  -  (x  -  Xl)(y  -  y2)  ’  * 

du  dv  dw 

+  7~™  =  o  ■  («) 


Here  x,  y,  and  z  represent  the  position  of  grid  points,  and  the  subscripted 
variables  (xj,  y*,  zt,  and  x2,  y2,  z2)  are  the  radars’  coordinates.  The  variance, 
(cr2),  of  a  dependent  variable  measures  the  expected  difference  between  a  sample 
and  the  true  value  of  that  variable.  The  error  variance  of  u,  v,  and  w  are  calculated 
at  each  grid  point.  In  Eqs.  Cl,  C2,  and  C3,  the  measurements  of  Vj,  V2t  and  W 
are  assumed  independent.  Thus,  in  calculating  error  variances,  the  covariances  can 
be  ignored  (Ray,  1980).  The  error  variances  for  u,  v,  w,  W,  Vj,  and  V2  are  given 


N  ,  ,2 

E(un  ~  U) 

N  -  1 


m 
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2  ^(Vn~v)2 

^  N-l  ’ 

n-l 

(C5) 

2  f  K-w)2 
-Z.  N  —  l  ’ 

n—l 

(CO) 

r  2  y>(^n-^)2 

°w  N- 1  ' 

n—l 

(Cl) 

2  y>(^ln-Vi)2 

*  S  iv-1  * 

(CO) 

2  ~  V2)2 

<’v’  N-l  ' 

n—l 

(C9) 

First  examine  the  u  component  equation.  Since  V\ ,  V2 ,  and  W  are  independent, 
only  the  derivation  of  the  first  term  need  be  examined  in  this  equation,  as  the 
mathematical  process  applied  will  give  the  corresponding  results  for  the  other  terms, 
as  well  as  the  v  component  equation.  The  first  term  of  the  u  component  is  given  in 
Eq.  Cl  as 


u  _ _ R\V\{y-yi) _ 

(i  -  xx)(y  -  y2)  -  (x  -  x2)(y  -  yt) 
Now  the  error  variance  for  this  term  is  found: 


(C 10) 


(uw  -  u)2  _  [fZiVjn(y  -  i/2)  -  RiVijy  -  t/2)]2 

N  ~ 1  ((a;-a?i)(y-y2)-(*-a?2)(y-yi)]2(iV-l) 

Equation  Cll  can  be  rewritten  in  the  form  of 


[Ri(y  ~  V2)? 


N  (Vxn-Vr)7 


1 


N—l 


[(x  -  xi)(y  -  y2)  -  (x  -  x2)(y  -  yi)]2 


(C 11) 


1 


(C 12) 
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which  is  the  same  as 


[Rt(  t/-y2)]W) 

l( x  -  xx )(y  -  y2)  -(x-  x2){y  -  yi)]2 


(C 13) 


This  is  the  expression  for  the  error  variance  for  the  first  term  of  u,  applying 
this  method  to  the  rest  of  Eq.  Cl  gives  the  error  variance  of  u  as 


_  (-Kip/  ~  +  \R*(v  ~  Vi)]  VO 

[(*  -  *i)(y  -  ya)  -  (*  -  *2)(y  -  yi)]2 


+  {[(•*  ~  *2 )(y  -  yi)  -  (*  -  zi )(y  -  y2)]2}(<7w2) 

1  [(x  -  Xx)(y  -  y2)  -  (x  -  x2)(y  -  yi)}2 

Correspondingly,  the  equation  for  the  error  variance  in  v  is  given  by 


(Cl  4) 


2  _  [fl2(a  ~  X2)]2((TV22)  +  [iZi(x  -  xx)]2{aVl2) 
[(x  -  x2)(y  -yx)-(x-xx )(y  -  y2)]2 


{[(*  ~ faXf  -xx)-(z-  zx)(x  -  x2))2}(aw2) 
[(x  -  x2)(y  -  yx)  -  (x  -  xx)(y  -  y2)]2 


(C 15) 


The  vertical  velocity  variance  is  found  by  taking  Eq.  C3  (the  anelastic  continu¬ 
ity  equation)  and  expressing  it  in  finite  difference  form.  In  this  form,  the  centered 
finite  differences  are  used  in  the  horizontal  (with  two  grid  spaces)  and  the  vertical 
(with  one  grid  space).  The  horizontal  velocities  are  at  level  k.  Since  we  are  in¬ 
tegrating  downward  the  value  of  to  at  the  top  of  the  layer  of  integration  tnust  be 
known  as  (k  +  ^).  Equation  C3  becomes 
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Appendix  D 


VAD  Analysis 


Assume  that  a  ground-based  radar  is  located  at  the  origin  of  a  Polar  coordinate  sys¬ 
tem.  The  kinematic  fields  from  the  Fourier  coefficients,  based  on  the  Velocity  Azimuth 
Display  (VAD)  technique  by  Browning  and  Wexler  (1968),  are  rederived  *using:  the  con¬ 
ventional  definition  of  the  azimuth  angle  (6)  in  a  clockwise  manner  (Browning  and  Wexler 
(1968)  used  counterclockwise);  the  convention  of  negative  Doppler  velocity  (Vj)  as  an  in¬ 
dication  of  flow  toward  the  radar;  and  the  fact  that  the  effects  of  the  beam  elevation  (<f>,  is 
a  small  angle)  are  negligible.  Assuming  a  uniform  wind  field,  Doppler  velocity  is  defined 


Vj  =  Vxcos^sin#  +  K,cos^cos  9 


(Dl) 


Here,  Vx  is  positive  towards  the  east  and  Vv  is  positive  towards  the  north.  The  velocity 
components  Vz  and  Vj,  can  be  expressed  in  terms  of  their  value  at  the  center  of  the  scanned 
circle  (subscript  0)  and  mean  linear  velocity  (overbarred  variables)  by 


1/  T/  .dV* 

K = v*° + it  +  - 

i,  ,  wv  w: 

V*  =  V*0  +  l£X  +  -d yV  • 


(D2) 


(D3) 


/ 
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Putting  Eqs.  D2  and  D3  into  Eq.  Dl,  using  x  =  rsin#  and  y  =  rcos#,  where  (r)  is  the 
range,  gives 


dV  dV 

Vm  =  Vzocosd)  +  rcos<f>-—^-sin20  -f  rcos<£-r-^cos#sin0 

ox  dy 


dVt 


dVu 


+  Vy0costf>  +  rcos<j>-~cos20  +  rcos<f>-^-cos0sm0 


dx 


dy 


(DA) 


Using  several  trigonometric  identities:  cos20  =  ^(1  +  cos20),  sin2#  =  *(1  —  cos2#), 
and  cos#sin#  =  ^sin2#,  Eq.  D4  becomes 


Vi 


r  f  dVz  dV  \ 

—  2COS<^\~g~  ~d^ )  ^*°cos<^s‘11^  +  Uyocos^cos^ 

dV, 


+  -cos<£ 


(dVz 

V  dy 


^sin2#  -  ^cos <f)(  —  ^t^cos2# 


+  St  )°in2e  -  2CUS\  dx  dy) 


(135) 


Equation  D5  can  be  put  into  a  Fourier  series  of  the  form 


1  °° 

Vd  =  -a0  +  ^2  (o„cosn0  -f  6„sinn#) 


(D  6) 


n=l 


This  gives  the  non-zero  Fourier  coefficients  of 


Vdvx  Wv\ 
00=^0,^—  +  —j. 

(Dl) 

°1  =  Uyocos^, 

(DS) 

6i  =  Uxocos^, 

(D9) 

/ 
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(DIO) 

(Dll) 


The  two-dimensional  horizontal  divergence,  stretching  deformation,  and  shearing  defor¬ 
mation  are,  respectively,  given  by  the  terms 


\  dx  dy  )  ’  \  dx  dy  dy  dx  ) 
Thus,  we  can  obtain  the  kinematic  properties  of  the  wind  from 


\V\  = 


(«;+*?)* 

COS  <t> 


(dV^.dVA 

\  dy  ^  dx  ) 


2  b2 
rcos<f> 


(D 12) 

(013) 

(014) 

(015) 


Neither  of  the  deformations  gives  individually  the  desired  total  rate  and  sense  of  the 
deformation  in  the  motion  unless  one  of  the  quantities  is  equal  to  zero.  To  derive  the  full 
value  from  this  concept,  the  two  types  of  deformation  axe  combined  to  give  the  resultant 
direction  of  elongation  (which  is  along  the  axis  of  dilatation)  and  the  net  rate  of  deformation 
(DEF)  (Saucier,  1965).  First,  the  angle  (a)  between  the  axis  of  dilatation  and  the  initial 
X  axis  is  calculated  using 
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av: ,  ,  , 

tan2a  =  ^i-+A  =  -2. 
a  14  014 


(me) 

_ «2 

a*  ay 

If  both  of  the  deformations  axe  of  the  same  sign,  the  axis  of  dilatation  is  rotated 


counterclockwise  through  the  angle  a.  Once  this  angle  has  been  determined  the  magnitude 
of  the  total  deformation  (DEF)  can  be  obtained  from 


DEF  = 


(£>17) 


Farther  explanation  of  these  terms  is  found  in  Chapter  3. 


